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EXECUTIVE  SUMMARY 


Four  Gaussian  puff  models  were  tested  with  seven  field  data  sets  collected  on  the  dispersion  of  fog 
oil  at  two  sites.  Three  data  sets  were  obtained  in  flat  terrain  at  Dugway  Proving  Ground  and  four  data  sets 
were  acquired  in  generally  flat  terrain  at  Camp  Atterbury  near  Columbus,  Indiana.  Data  on  average 
concentration  over  the  time  of  release  were  obtained  at  distances  of  25  m  -  200  m  at  Dugway  and  50  m  - 
675  m  at  Camp  Atterbury.  Data  were  acquired  over  several  transects  across  the  plume  downwind  from  the 
source.  In  support  of  the  in-plume  average  concentration  measurements,  ihere  were  time-dependent 
source  measurements  made  along  with  time-dependent  vertical  profiles  of  meteorological  data  (mean 
values  and  turbulence  quantities). 

The  models  tested  were  the  BEAR,  INPUFF-Onsite,  INPUFF-PG,  and  RIMPUFF.  These  models 
treat  the  release  as  a  series  of  Gaussian  puffs  which  get  transported  and  dispersed  downwind  with  a  time- 
dependent  wind  field.  The  results  show  that  the  four  Gaussian  puff  models  tested  can  predict  within  a 
factor  of  2-3  under  the  convective  conditions  tested  to  distances  of  about  250  m.  Beyond  that  distance, 
the  plume  tends  to  rise  leading  the  models  to  significantly  overpredict  average  concentrations  at  ground 
level.  No  data  were  acquired  at  these  two  sites  under  stable  conditions.  (Fog-oil  plumes  were  measured 
under  stable  conditions  at  the  Meadowbrook  site  in  Red  Bluff,  California  but  model  testing  with  those  data 
is  not  the  subject  of  this  report. ) 

One  difficulty  with  the  application  of  the  above  four  models  to  the  seven  data  sets  is  the  simple 
nature  of  the  Gaussian  models  tested.  Indeed,  simple  models  are  needed  for  the  practical  and  rather 
routine  smoke  plume  modeling  used  by  the  U.S.  Army.  However,  simple  models  usually  require  simple 
inputs  such  as  a  single-level  wind  speed  as  a  function  of  time.  Should  the  wind  speed  at  the  release 
height  be  used  or  is  the  10-m  value  more  representative?  An  entire  profile  of  such  meteorological  data  is 
available  from  the  experiments.  What  is  the  initial  puff  radius  and  height  above  the  ground  once  the  initial 
momentum  of  the  plume  exiting  the  three  release  ports  is  accounted  for?  There  are  no  definitive  answers 
to  such  questions  (known  at  this  iime)  and  the  models  are  sensitive  to  the  choices  made.  In  any  case,  the 
agreement  between  tne  modeis  and  data  is  reasonable  for  such  flat  terrain  applications  and  distances 
equal  to  or  less  than  about  250  m. 

Although  the  Gaussian  models  tested  with  the  fog-oil  data  indicated  an  overall  agreement  with  the 
data  within  a  factor  of  2-3  for  the  average  concentrations  measured,  there  were  problems  with  the  models 
(under  these  largely  unstable  atmospheric  conditions)  in  that 

(i)  such  agreement  was  only  within  the  first  250  m  from  the  smoke  generator  with  more  significant 
discrepancies  at  further  distances, 

(ii)  there  is  a  tendency  to  significantly  overpredict  or  underpredict  average  concentrations  at 
distances  less  than  100  m  and  to  significantly  overpredict  at  distances  greater  than  about  250  m.  The 
discrepancies  found  at  longer  distances  (greater  than  about  250  m)  is  due  to  a  rise  in  the  plume 
centerline,  and 

(iii)  the  decay  of  concentration  with  distance  as  predicted  by  the  models  does  not  agree  with  the 
data  and  is  likely  to  be  the  result  of  the  rising  centerline  seen  in  the  observed  plume.  The  rising  centerline 
phenomenon  has  been  observed  in  other  field  and  laboratory  studies  and  cannot  be  predicted  by  these 
Gaussian  models. 

In  a  companion  report,  a  stochastic  model  is  compared  to  the  Camp  Atterbury  data  and  is  shown  to 
predict  the  rise  of  the  plume  under  convective  conditions  quite  well  along  with  the  correct  decay  of 
concentration  with  distance.  The  lack  of  treatment  of  convective  turbulence  in  the  Gaussian  puff  models 
is  apparently  the  cause  of  the  discrepancies  between  the  model  predictions  and  the  data  beyond  about 
250  m.  Future  work  will  involve  the  testing  of  complex  terrain  versions  of  these  types  of  models  with  both 
the  stable  and  unstable  plume  dispersion  data  at  the  Meadowbrook  site. 
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1.0  INTRODUCTION 


This  report  provides  an  evaluation  of  four  of  the  more  promising  mathematical  models  for  the 
prediction  of  the  dispersion  of  fog-oil  smoke  discharged  from  M3A3E3  or  M3A4  smoke  generators.  The 
data  base  used  for  model  testing  consists  of  seven  data  sets:  three  field  trials  conducted  at  Dugway 
Proving  Ground,  Utah  during  March/April  1985  (Liljegren  et  al.,  1988)  and  four  field  trials  conducted  at 
Camp  Atterbury,  Indiana  in  November  1987  (Liljegren  et  al.,  1989). 

The  purpose  of  the  model  evaluation  portion  of  the  current  research  program  is  to  identify  the 
superior  modeling  approaches  and  individual  models  for  fog-oil  smoke  dispersion.  Some  model 
improvement  may  be  necessary  to  satisfy  the  objectives  of  the  U.S.  Army.  Such  validated  model(s)  could 
then  be  used  in  two  ways: 

(a)  In  environmental  impact  analyses  to  determine  the  health  and  environmental  effects  of  smokes 
used  in  training  exercises  carried  out  by  the  U.S.  Army.  Of  greatest  interest  here  is  the  prediction  of 
dosage  (time-integrated  concentration),  particle-size  distributions,  and  deposition  rates  at  distances  just 
downwind  of  the  generator  (e.g.,  20m)  to  several  kilometers  downwind,  and 

(b)  During  the  execution  of  training  exercises  with  smokes.  In  such  circumstances,  the  model(s) 
would  be  used  in  a  real-time  mode  to  determine  whether  a  specific  fog-oil  trial  should  be  undertaken  given 
current  meteorological  conditions.  The  prediction  of  the  extent  of  the  visible  plume  would  be  useful  as 
well  as  the  quantities  identified  in  (a). 

The  measurement  program  undertaken  at  Dugway  Proving  Ground  and  Camp  Atterbury  was  aimed 
at  providing  the  data  in  (a).  Color  video  photographs  of  the  plume  were  taken  at  both  sites;  however, 
plume  visibility  was  not  an  objective  of  either  set  of  field  tests  and  so  only  qualitative  information  on  plume 
visibility  was  obtained. 

An  early  evaluation  of  fog-oil  plume  models  was  carried  out  as  part  of  the  current  program  (Policastro 
et.  al.,  1985).  In  that  work,  four  models  were  tested  with  three  data  sets  from  Smoke  Weeks  III  and  IV.  The 
models  were  Act  II  (Sutherland  et  al.,  1982  and  Petraska,  1984),  MAD  PUFF  (Petraska,  1984),  COMBIC 
(Hoock  et  al.,  1982),  and  Ludwig  (Ludwig,  1977).  Although  there  was  not  a  large  difference  among  the 
predictions  of  the  models,  it  became  clear  that  models  that  could  handle  time-varying  meteorology  would 
be  the  most  useful  considering  the  dual  purposes  (see  (a)  and  (b)  above)  that  the  modeling  must  satisfy. 
The  model/data  comparisons  provided  in  that  paper  (Policastro  et.  al.,  1985,  and  also  Appendix  A  of  this 
report)  could  not  lead  to  any  definitive  conclusions  because: 

(i)  The  Smoke  Week  III  and  IV  data  for  dosage  were  adjusted  subjectively  by  Dugway  personnel  in 
order  to  provide  more  physically  reasonable  data.  The  amount  of  this  correction  or  adjustment  is  unknown 
and  casts  doubt  as  to  the  accuracy  of  that  data  base,  and 

(ii)  The  data  were  acquired  at  distances  only  within  about  90  m  from  the  source  and  involved 
generator  run  times  of  only  a  few  minutes. 

The  three  data  sets  acquired  during  the  March/April  1985  time  period  at  Dugway  Proving  Ground 
and  the  four  data  sets  acquired  during  the  November  1987  time  period  at  Camp  Atterbury,  Indiana  provide 
new  sources  of  data  for  model  testing.  Measurements  were  made  to  slightly  longer  distances  (out  to  200 
meters  at  Dugway  and  out  to  675  meters  at  Camp  Atterbury)  and  for  generator  run  times  close  to  one 
hour.  No  adjustments  were  made  to  these  data  in  any  way. 

Presented  below  is  a  description  of  the  three  models  tested  with  these  data  and  the  results  of  the 
model/data  comparisons. 


2.0  THE  MODELS 

An  examination  of  the  literature  reveals  that  there  are  three  general  modeling  approaches  that  can  be 
used  to  simulate  the  dispersion  of  fog-oil  plumes.  These  are: 
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(a)  Gaussian  puff  method 


In  this  treatment,  a  continuous  release  is  divided  up  into  a  number  of  puffs  whereby  each  puff  is 
advected  and  dispersed  downwind  based  on  time-dependent  meteorology.  Models  within  this  category 
differ  based  largely  upon  the  dispersion  coefficients  used  as  well  as  wind  speeds  (wind  field)  employed. 
This  technique  has  the  advantage  of  handling  a  variable  source  release  rate  which  is  common  with  U  S. 
Army  smoke  generators.  Models  in  this  category  are  Ludwig  (Ludwig,  1977),  RIMPUFF  (Mikkelsen  et  al. , 
1984),  and  INPUFF  (Petersen,  et  al.  1984).  Note,  the  Ludwig  model  is  also  known  as  the  BEAR  model. 
This  approach  is  amenable  to  handling  time-dependent  meteorological  data  from  a  tower  at  the  site  of 
interest. 

(b)  Gaussian  plume  method 

This  methodology  assumes  that  the  continuous  Gaussian  plume  formulation  can  best  represent  the 
fog  oil  release.  The  treatment  of  long  release  times  is  handled  through  the  dependence  of  the  dispersion 
coefficients  (ax,  ay  and  a2)  with  release  time.  For  large  release  times,  for  instance,  the  ay  coefficient  is 
large  leading  to  a  wider  plume.  That  same  plume  might  just  as  well  have  been  predicted  by  the  Gaussian 
puff  method  with  many  puffs  being  spread  out  in  different  wind  directions  over  that  long  release  time,  with 
each  puff  having  a  small  ay.  It  is  not  a  priori  clear  whether  the  Gaussian  plume  method  can  be  as  accurate 
as  the  Gaussian  puff  method  for  the  long  release  times  expected  for  fog-oil  plumes.  Models  that  are 
represented  in  the  category  of  Gaussian  plume  method  are:  ACT  II  (Sutherland  et  al.,  1982  and  Petraska, 
1984),  MAD  Puff  (Petraska,  1984),  COMBIC  (Hoock  et  al.,  1982),  MSMOKE  (Hansen,  1984),  HZRD 
ll(Pennsyle,  1984),  and  HECSMOKE-I  (Cheney  and  Dumbauld,  1979). 

(c)  Monte-Carlo  method 

This  methodology  represents  the  fog-oil  release  as  the  emission  of  large  numbers  of  particles  that 
are  advected  and  dispersed  with  the  wind.  The  advantage  to  this  method  is  that  it  is  better  suited  to 
handling  the  more  fundamental  aspects  of  atmospheric  turbulence  (does  not  rely  on  Gaussian  dispersion 
coefficients  as  do  the  other  two  methods)  and  can  more  accurately  treat  surface  boundary  conditions.  This 
method  is  more  suitable  in  complex  terrain  situations  as  well.  Unfortunately,  the  existing  models  in  this 
category  must  still  be  considered  as  research  tools  and  are  only  now  entering  the  competitive  arena  with 
the  Gaussian  puff  and  plume  models.  Models  in  this  category  are  MoCaPD  (Huang  and  Frost,  1982), 
Dunn-Boughton  (Boughton  and  Dunn,  1983)  and  Schorling  (Schorling,  1986). 

2.1  Models  Evaluated 

For  the  purposes  of  this  evaluation  of  models  with  the  Dugway  and  Camp  Atterbury  data,  three 
Gaussian  puff  models  have  been  chosen  for  testing.  They  are  BEAR,  INPUFF,  and  RIMPUFF.  These 
models  satisfy  the  criteria  of  handling  time-dependent  meteorological  data  and  ease  of  use. 

The  simplest  of  the  three  models  is  the  BEAR  model.  This  model  permits  as  input  time-dependent 
meteorological  data  for  use  in  dispersing  puffs.  Each  time  that  meteorological  inputs  are  read,  a  new  time 
step  for  the  release  of  puffs  is  determined.  The  model  permits  the  simulation  of  changing  wind  speeds 
and  directions  during  the  period  of  smoke  generation.  As  a  result,  plume  meander  is  treated  by  means  of 
this  time-dependent  formulation.  The  Gaussian  puff  dispersion  coefficients  employed  in  the  model  are 
the  Pasquill-Turner  coefficients  that  are  dependent  on  stability  class.  Although  developed  for  continuous 
plumes,  they  are  applied  by  Ludwig  to  puffs.  Meteorological  data  are  read  into  the  model  in  three  minute 
time  steps,  consistent  with  the  averaging  time  of  the  Pasquill-Gifford-Turner  dispersion  coefficients 
(Turner,  1970). 

The  RIMPUFF  model  is  similar  in  its  formulation  to  the  BEAR  model.  However,  the  dispersion 
coefficients  are  functions  of  downwind  distance  based  on  either: 

(a)  the  Pasquill-Gifford-Turner  stability  class,  or 

(b)  the  variances  in  the  wind  direction  using  the  work  done  by  Smith  and  Hay  (1961). 


In  this  model,  meteorological  time  steps  are  typically  15  seconds  with  the  dispersion  coefficients 
adjusted  to  that  time  scale.  In  this  way,  plume  meander  is  simulated  directly  from  the  variation  of  wind 
direction  with  time.  Meander  is  not  included  in  the  dispersion  coefficients. 

The  INPUFF  model  is  also  similar  to  BEAR  and  RIMPUFF  in  that  it  divides  the  release  into  a  series  of 
puffs.  For  smoke  applications,  the  meteorological  time  step  is  typically  2-3  minutes.  The  puff  release  rate  is 
independently  set  and  may  be  chosen  to  be  once  every  10  seconds.  This  dichotomy  parallels  that  of  the 
Ludwig  model.  Three  dispersion  algorithms  are  used  within  INPUFF  for  dispersion  downwind  of  the 
source.  The  user  may  use  the  Pasquill-Gifford-Turner  (P-G)  scheme  or  the  On-Site  scheme,  so  named 
because  it  requires  specification  of  the  variances  of  the  vertical  and  lateral  wind  direction.  The  On-Site 
scheme  is  a  synthesis  of  work  performed  by  Draxler  (Draxler,  1976)  and  Cramer  (Cramer,  1976).  The  third 
dispersion  scheme  is  for  long  travel  times  in  which  the  growth  of  the  puff  becomes  proportional  to  the 
square  root  of  time. 

The  treatment  of  the  wind  field  for  a  typical  smoke  release  is  handled  in  the  following  way  with  these 
models.  The  BEAR  model  requires  the  input  of  the  10-m  level  wind  speed,  wind  direction,  and  relative 
humidity  as  a  function  of  time.  The  model  handles  only  one  time  series  of  measurements  at  only  one  level. 
The  BEAR  model  requires  the  meteorological  data  to  be  input  as  roughly  three-minute  averages,  in 
concert  with  the  averaging  time  of  the  dispersion  coefficients. 

For  RIMPUFF,  meteorological  input  is  more  general.  That  model  can  handle  multiple  meteorologicr 
stations;  however,  each  station  must  have  its  time  series  of  data  provided  at  only  one  level.  The  data, 
however,  may  be  taken  at  different  reference  heights  whereby  the  model  interpolates  all  to  a  single 
inputted  reference  height.  The  input  height  of  the  meteorological  measurements  for  this  model  is  not 
fixed  at  10  meters  as  in  the  other  models,  rather  any  height  can  be  specified  and  the  model  will  create  a 
vertical  wind  profile  from  the  chosen  reference  height.  The  model  computes  the  u  and  v  wind 
components  on  a  grid;  1/r2  interpolation  in  space  and  linear  interpolation  in  time  is  used  to  provide  the 
wind  speed  components  at  each  grid  point  at  each  meteorological  time  step.  The  wind  field  created  in  this 
manner  will  not  necessarily  be  divergence  tree. 


For  INPUFF,  the  user  must  supply  the  wind  speed  and  direction  within  a  grid  system.  The  user 
defines  the  coordinate  and  size  of  each  grid  square  as  well  as  the  extent  of  the  meteorological  region.  The 
meteorological  region  may  be  larger  or  smaller  than  the  modeling  region.  The  meteorological  data  (wind 
speed,  wind  direction,  Oq,  Cty  ,  mixing  height,  and  stability  class)  are  read  in  as  a  function  of  time.  The  P-G 
scheme  does  not  use  Gq  and  G<j,  values,  while  the  On-Site  scheme  does  not  use  the  stability  class  A 
meteorological  preprocessor  definea  by  the  user  may  compute  (single-layer)  winds  and  the  output  of  that 
preprocessor  may  be  input  into  INPUFF. 

For  a  single  meteorological  tower  supporting  a  smoke  experiment,  all  three  models  can  have  their 
wind  field  treatment  reduce  to  an  identical  formulation;  the  10-m  wind  level  as  a  function  of  time.  The 
averaging  period  for  that  10-m  level  data  will,  of  course,  be  different  depending  upon  the  model  due  to 
the  needs  of  its  dispersion  scheme. 

For  the  Dugway  and  Camp  Atterbury  tests,  the  models  being  tested  are  similar  in  many  ways  with  the 
main  exception  being  in  their  treatment  of  the  dispersion  of  the  puffs.  For  each  of  the  models,  coding 
changes  were  made  to  predict  dosage  at  each  receptor.  Each  model  had  predicted  concentration  with 
time.  Table  2.1  gives  a  comparison  of  the  theoretical  features  of  these  three  models. 

3.0  THE  DUGWAY  FIELD  DATA 


This  section  provides  a  short  summary  of  the  data  acquired  in  the  March/ April  1985  field  tests  at  Dugway 
Proving  Ground  (Liljegren  et  al.,  1988).  Three  data  sets  are  used  for  model  testing  in  the  next  section. 
More  details  concerning  all  the  data  acquired  as  well  as  the  equipment  and  methodology  used  is  given  in 
Liljegren  et  al.(1988). 
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Table  2.1  Comparison  of  theoretical  features  of  the  INPUFF,  BEAR,  and 
RIMPUFF  Models 


INPUFF 

REAR 

RIMPUFF 

Wind  Field 

provided  by  user, 
variable  in  space 
and  time 

single  point  and 
elevation  but 
variable  with  time 

multiple  station 
measurements  in¬ 
terpolated  >n 
space  and  time 

Dispersion 

Methodology 

Gaussian  puff 

Gaussian  puff 

Gaussian  puff 

Dispersion 

PGT  stability 

PGT  stability 

PGT  stability 

Coefficients 

classes  (f)  or  on-site 
scheme  (a) 

classes(f) 

classes(f)  or  on-site 
scheme  (b) 

Meteorological 
Time-Averaging 
Period  (c) 

2-5  minutes 

3  minutes 

15  seconds 

Puff  Release 

Interval  (d) 

1  second 

1  second 

1  second 

* 

Treatment  of 

Wind  Shear  on 

Puff  Dispersion  (e) 

No 

No 

Yes 

Source  Conditions: 

Time  Dependent 

Yes 

No 

Yes 

Multiple  Sources 

Yes 

No 

Yes 

(a)  The  on-site  scheme  uses  Oe  and  measurements  in  formulas  for  Oy  and  az  derived  by  Irwin 
(1 983)  from  a  combination  of  the  Draxler  (1976)  and  Cramer  (1976)  formulas. 

(b)  The  on-site  scheme  here  employs  a9  numbers  in  the  Smith-Hay  dispersion  formulas  (1961). 

(c)  The  times  provided  here  are  recommended  by  the  modelers  based  on  the  dispersion  coefficients 
used  in  the  model.  The  INPUFF  (PG  scheme)  and  the  BEAR  Model  use  the  same  dispersion 
coefficients. 

(d)  The  one  second  values  are  the  lowest  values  that  are  recommended  by  the  modelers.  Larger 
values  reduce  computer  time  but  also  reduce  prediction  accuracy  dependent  upon  the  extent  of 
departure  from  the  1  second  values. 

(e)  Ignoring  the  effect  of  wind  shear  on  individual  puff  dispersion  leads  to  a  Gaussian  distribution  in 
the  vertical  direction.  The  RIMPUFF  Model  does  not  yield  a  Gaussian  distribution  in  the  vertical 
direction  as  a  result  of  its  treatment  of  vertical  wind  shear  on  plume  dispersion. 

(f)  PGT  stability  classes  refer  to  Pasquill-Gifford-Turner  stability  classes  A-F,  where  A  represents  the 
most  unstable  conditions  and  F  the  most  stable  conditions. 
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Figure  3,1  M3A3E3  Smoke  Generator  assemblies  and  major  components 
( from  M3A3E3  Operations  Manual). 
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nd  release  rate  (g/s)  as  a  function  of  time  for  test 


Figure  3.3.  Location  of  test  site  with  respect  to  Horizontal  and  West  Vertical  grids  at 
Dugway  Proving  Ground.  The  sampling  network  is  located  in  the  shaded 
portion  of  the  figure. 


All  three  tests  employed  a  M3A3E3  military  smoke  generator  (see  Figure  3.1).  It  produces  smoke  by 
vaporizing  fog  oil  and  ejecting  it  into  the  atmosphere.  Time-dependent  measurements  of  exit  temperature 
were  taken  by  instrumenting  each  of  the  three  exit  nozzles  on  the  generator  with  a  chromel/alumel  (type  k) 
thermocouple.  Estimates  of  the  exit  velocity  of  the  smoke  were  obtained  periodically  during  a  test  using  a 
pitot  tube.  The  calculation  of  the  exit  velocity  accounted  for  the  high  temperature  of  the  exhaust  stream  as 
well  as  the  increased  density  of  the  exhaust  stream  due  to  the  fog  oil.  A  scale  was  used  to  provide  the 
remaining  weight  of  fog  oil  in  a  55  gallon  drum  as  a  function  of  time. 

A  summary  of  the  data  acquired  at  the  source  for  the  three  experiments  is  given  in  Table  3.1  with  an 
illustration  of  the  time-dependent  source  data  given  in  Figure  3.2.  Note  that  the  release  rate  varied 
significantly  during  the  test.  This  highlights  the  need  for  time-dependent  mass  release  data  and  a  model 
that  can  treat  that  time-dependent  release. 

The  meteorological  data  were  collected  from  two  32-m  towers  located  at  the  Dugway  test  site  (  see 
Figure  3.3  ).  The  test  site  for  the  fog-oil  smoke  release  was  located  between  the  meteorological 
instrument  towers  on  Horizontal  Grid  and  West  Vertical  Grid.  The  site  of  the  testing  passes  the 
micrometeorological  test  of  uniform  terrain  in  all  directions  to  a  distance  equal  to  100  times  the  tower 
heights.  The  roughness  length  for  the  site  is  2  cm.  The  primary  vegetation  for  the  site  is  a  desert  shrub 
having  a  height  of  7  to  30  cm  arid  spaced  at  one-half  to  one  meter  intervals. 

At  each  tower,  wind  speeds  were  measured  with  cup  anemometers  at  2,  4,  8,  16,  and  32m. 
Fluctuations  in  the  horizontal  and  vertical  wind  angles  were  measured  with  bivanes  mounted  at  4,  8,  16, 
and  32m.  Horizontal  wind  fluctuations  at  the  2m  level  were  measured  with  a  direction  vane.  On  both 
towers,  temperatures  were  measured  at  the  2,  4,  8,  16,  and  32m  levels.  The  meteorological  data  on  both 
towers  were  available  at  one  second  intervals.  A  summary  of  the  meteorological  data  acquired  during  the 
three  tests  is  given  in  Table  3.2.  Considering  that  only  three  test  data  sets  were  to  be  compared  with 
model  predictions,  a  comparison  of  the  acquired  meteorological  data  with  long-term  stability  class  /  wind 
speed  data  at  Dugway  Proving  Ground  was  not  found  to  be  fruitful. 

The  location  of  the  sampling  network  within  the  West  Vertical  Grid  is  illustrated  by  the  shaded  region 
in  Figure  3.4.  A  high  resolution  network  was  laid  out  in  the  shaded  area  of  Figure  4  between  rows  2  and  6 
and  columns  J  and  N  (  see  Figures  3.4  and  3  5).  As  illustrated  in  Figure  3.5,  this  network  is  essentially 
diamond-shaped.  It  is  composed  of  transects  at  25  m,  50  m,  100  m,  200  m,  and  400  m  and  is  symmetric 
about  the  200  m  transect.  Eleven  samplers  were  located  on  both  the  25  m  and  50  m  transects  at  spacings 
of  5  m  and  10  m  respectively.  Nine  samplers  were  located  on  the  100  m,  200  m,  and  sampling  plan,  the 
smoke  source  could  be  located  near  J4(0,0)  or  N4(400,0)  depending  on  the  wind  direction.  At  each 
sampler  placed  on  operation,  total  mass  of  fog  oil  collected  was  measured.  More  details  about  the 
sampling  system  appear  in  the  Dugway  data  report  (Liljegren  et  al,  1988). 

Average  concentration  for  test  cases  T0090,  T0010,  and  T0011  are  summarized  in  the  isopleth 
plots  presented  in  Figures  3. 6, 3. 7,  and  3.8.  Note  that  for  test  T0011  (11  April  1985),  the  wind  direction 
changed  while  the  generators  were  being  field  repaired  after  breakdown.  To  accommodate  this  wind  shift, 
the  generator  was  relocated  prior  to  restarting  the  test;  both  locations  are  indicated  in  Figure  3.8.  Model 
predictions  will  be  compared  to  the  data  on  a  transect  basis  in  the  next  section  of  this  report. 

4.0  COMPARISON  OF  MODELS  WITH  FIELD  DATA  TAKEN  AT  DUGWAY 

PROVING  GROUND 

Application  of  the  BEAR,  RIMPUFF,  and  INPUFF  models  to  these  three  data  sets  was  accomplished 
in  a  straight-forward  manner.  The  meteorological  data  for  these  three  data  sets  were  available  on  a  one 
second  basis  from  Dugway  Proving  Ground  from  which  time  averages  were  prepared  based  on  the 
appropriate  averaging  of  meteorological  data  required  by  each  model.  Next,  the  time  series  of  mass 
release  of  fog  oil  was  input  for  both  schemes  of  the  INPUFF  model,  with  a  time  averaging  period  of  3 
minutes.  The  RIMPUFF  and  BEAR  models  employ  only  one  (average)  mass  release  rate.  Both  the  BEAR 
and  the  INPUFF  models  used  a  time  averaging  period  of  3  minutes  for  the  meteorological  data,  while  the 
RIMPUFF  model  used  a  time  averaging  period  of  15  seconds. 
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Table  3.2.  Summary  of  meteorological  data  for  the  three  Dugway  data  sets 


Notes:  1  All  values  are  averages  over  th  a  period  of  the  test. 

2.  All  values  are  for  the  10m  level  except  for  the  relative 
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Figure  3.4.  Layout  of  main  sampling  network  and  coordinate  sytem.  Filled  circles 

indicate  sampler  locations.  Diamonds  (0)  indicate  sources.  Shaded  areas 
locates  High  Resolution  Network  illustrated  in  Figure  3.3 
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Figure  3.5  Higti  resolution  sampling  network  indicated  by  the  shaded  area  in  Figure 
3  3.  t  illod  circles  maiked  the  locations  of  concentration  samplers 


Figure  3.6.  Isopleths  of  average  concentration  for  test  T0009.  Sampler  locations  are 

indicated  by  small  filled  circles.  The  arrow  at  the  source  location  indicates  the 
mean  wind  direction;  its  length  is  proportional  to  the  mean  wind  speed.  The 
shaded  sector  around  the  arrow  indicates  the  mean  wind  direction;  its  length 
is  proportional  to  the  mean  wind  speed.  The  shaded  sector  around  the  arrow 
indicates  the  standard  deviation  of  the  wind  direction. 
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Figure  3.7.  Isopleths  of  average  concentration  for  test  T001 0.  Sampler  locations  are 

indicated  by  small  filled  circles.  The  arrow  at  the  source  location  indicates  the 
mean  wind  direction;  its  length  is  proportional  to  the  mean  wind  speed.  The 
shaded  sector  around  the  arrow  indicates  the  standard  deviation  of  the  wind 
direction. 
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Figure  3.8. 
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Isopleths  of  average  concentration  for  test  T001 1 .  Sampler  locations  are 
indicated  by  small  filled  circles.  The  arrow  at  the  source  location  indicates  the 
mean  wind  direction;  its  length  is  proportional  to  the  mean  wind  speed.  The 
shaded  sector  around  the  arrow  indicates  the  standard  deviation  of  the  wind 
direction. 
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Figure  4.1.  Comparison  of  RIMPUFF,  BEAR,  INPUFF-Onsite  and  INPUFF-PG 
mathematical  models  with  average  concentrations  from  test  T0009 
Dugway  Proving  Ground. 
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Figure  4.2.  Isopleths  of  average  concentration  for  test  T0009  at  Dugway  Proving 

Ground.  Sampler  locations  are  indicated  by  small  filled  circles.  The  arrow 
at  the  source  location  indicates  the  mean  wind  direction;  its  length  is 
proportional  to  the  mean  wind  speed.  The  shaded  sector  around  the  arrow 
indicates  the  standard  deviation  of  the  wind  direction. 
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Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  BEAR  model  for  test  T0009  at  Dugway  Proving  Ground. 


60001 


Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  INPUFF-Onsite  model  for  Test  T0009  at  Dugway  Proving  Ground. 
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gure  4.5.  Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  INPUFF-PG  model  for  Test  T0009  at  Dugway  Proving  Ground. 
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Figure  4.7.  Comparison  of  RIMPUFF,  BEAR,  INPUFF-Onsife  and  INPUFF-PG 
mathematical  models  with  average  concentrations  from  test  TOOIC 
Dugway  Proving  Ground. 
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RESTART: 


Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  BEAFl  model  for  test  T0010  at  Dugway  Proving  Ground. 


Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  INPUFF-Onsite  model  for  Test  T0010  at  Dugway  Proving  Ground. 
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Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  RIMPUFF  model  for  Test  T0010  at  Dugway  Proving  Ground. 
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Figure  4. 1 4.  Isopleths  of  average  concentration  for  test  T001 1  at  Dugway  Proving 
Ground.  Sampler  locations  are  indicated  by  small  filled  circles.  The 
arrow  at  the  source  location  indicates  the  mean  wind  direction;  its  length 
is  proportional  to  the  mean  wind  speed.  The  shaded  sector  around  the 
arrow  indicates  the  standard  deviation  of  the  wind  direction. 
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SMOKE:  FOG  OIL 

BEGIN:  9:38° 

STOP:  9:45 


Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  BEAR  model  for  test  T001 1  at  Dugway  Proving  Ground. 


Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  INPUFF-PG  model  for  Test  T001 1  at  Dugway  Proving  Ground. 
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Figure  4.18.  Isopleths  of  average  concentration  computed  on  a  rectangular  grid  using 
the  RIMPUFF  model  for  Test  T001 1  at  Dugway  Proving  Ground. 


For  INPUFF,  two  modeling  schemes  were  tested;  the  On-Site  scheme  and  the  P-G  scheme.  The 
on-site  scheme  employs  measured  turbulence  data  (Gq  and  a^)  whereas  the  P-G  scheme  requires  the 
input  of  a  Pasquill  stability  class. 

All  three  field  tests  were  taken  under  convective  conditions  with  the  exception  that  neutral 
atmospheric  conditions  occur  for  the  first  part  of  test  T001 1  a.  Color  video  photographs  recorded  a  rise  of 
the  plume  centerline  with  distance  beyond  about  300-500  m  downwind. 

The  complete  set  of  model/data  comparisons  are  presented  in  Figures  4.1  to  4.18.  Figures  4.1  to 

4.6  refer  to  model/data  comparisons  for  test  T0009.  Figures  4.7  to  4.12  refer  to  model/data  comparisons 

for  test  T001C.  Figure  4.13  to  4.18  refer  to  tesi  T0011.  It.  those  fiyures,  two  types  of  evaluations 

are  made: 

(a)  comparisons  of  average  concentrations  in  mg/m3  between  predicted  and  obser/ed  over  the  time 
of  the  smoke  release.  Presentations  of  these  pointwise  model/data  comparisons  are  made  by  transect 
downwind  of  the  source  for  each  experiment,  and 

(b)  comparisons  of  isocontours  of  dosage  between  data  and  model  predictions.  Contours  are 
typically  prepared  for  4,  8,  16,  and  32  mg/m3-  The  purpose  of  these  contour  comparisons  is  to  determine 
if  the  spatial  pattern  of  the  predicted  average  smoke  concentrations  compares  well  with  the  observed 
pattern. 

The  type  (a)  comparisons  represent  only  point-by-point  evaluations  of  model  predictions.  In 
actuality,  the  presence  of  wind-direction  errors  in  the  data  or  the  model  predictions  could  lead  to 
inaccurate  pointwise  comparisons,  yet  the  spatial  pattern  of  the  predictions  and  observations  may  look 
alike.  For  example,  the  spatial  patterns  may  look  alike  but  be  rotated  from  each  other.  In  such  a  case,  the 
pointwise  comparisons  may  look  poor  but  the  spatial  patterns  may  look  similar.  In  this  way,  the  type  (a)  and 
type  (b)  comparisons  complement  each  other  by  evaluating  the  models  point-by-point  and  through  their 
isocontour  patterns. 

For  the  input  meteorological  data,  the  10  m  height  was  used  to  represent  the  wind  data  required  by 
the  models.  This  10  m  level  for  the  meteorological  data  is  recommended  by  the  model  developers.  The 
smoke  release  height  is  at  2.0  m  from  the  ground;  therefore,  the  center  of  each  puff  remains  at  the 
elevation  of  release  (2.0  m)  during  its  entire  downwind  dispersion  history. 

It  was  observed  in  the  field  at  Dugway  Proving  Ground  that  the  plume  centerline  rises  under 
convective  conditions.  That  observation  contradicts  the  models'  assumption  that  the  centerline  of  the 
plume  remains  at  the  elevation  of  release  (for  a  nonbuoyant  source  like  the  fog  oil  release).  However,  the 
fog  oil  smoke  data  were  acquired  at  distances  close  enough  to  the  generator  (within  200m )  that  the  rising 
centerline  phenomenon  had  not  yet  had  a  chance  to  manifest  itself  in  these  data  sets.  As  will  be  seen  later 
in  the  Camp  Atterbury  comparisons,  the  rising  centerline  was  observed  in  the  data  for  distances  of  about 
250  -  400  m  and  beyond  under  the  convective  conditions  present  during  that  study. 

4.1  Discussion  of  the  Field  Data  for  the  Dugway  Tests 

An  examination  of  the  average  concentration  isocontours  (Figures  4.1,  4.7,  and  4.13)  reveals  a 
general  consistency  among  the  data  presented  for  tests  T0009,  T0010,  and  T0011.  The  three 
experiments  are  similar  in  terms  of  generator  run  time,  downwind  sampler  coverage,  and  meteorological 
conditions.  As  a  result,  the  dosages  measured  should  be  similar.  Indeed,  an  examination  of  Figures  4.1, 

4.7  and  4.13  reveal  that  contours  of  average  concentrations  of  4  to  64  mg/m3  are  present  in  each  test  for 
samplers  out  to  200  m  from  the  source. 

Other  areas  of  data  consistency  are  as  follows: 

(a)  The  lateral  distribution  of  average  concentrations  appears  to  be  Gaussian  in  shape  based  on  an 
examination  of  crosswind  measurements  for  each  of  the  transects.  The  result  is  consistent  with  the 
Gaussian  predictions  of  the  models. 
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(b)  Average  concentrations  decrease  in  a  systematic  way  in  the  longitudinal  direction.  The 
measured  values  appear  to  decline  by  factors  of  about  30-90  from  the  25  m  transect  to  the  200  m  transect. 
The  model  predictions  decline  about  the  same  amount,  roughly  by  factors  of  25-70.  The  decrease  in 
concentrations  with  distance  tends  to  be  in  a  geometric  manner  with  increasing  rates  of  decline  with 
distance.  For  example,  the  centerline  distance  to  8  mg/m3  should  be  greater  than  twice  the  centerline 
distance  to  the  4  mg/m3  average  concentration  value.  Since  in  dispersion  theory,  ay  and  Gz  inc.easpc:  at  a 
rate  faster  than  linear  as  a  function  of  x  for  stability  classes  B  and  C,  then  the  resulting  concentration 
computed  should  decrease  more  rapidly  than  linearly.  The  rate  of  decrease  in  both  models  and  data  is 
faster  than  linear  and  generally  of  the  same  order  of  magnitude  for  models  and  data. 

(c)  A  rtuticSdbie  meandering  of  tne  centeriine  ot  the  isopleths  is  observed  in  tests  T0009  and 
T0011.  No  such  meandering  seems  to  have  occurred  in  test  T0010.  The  distortion  is  seen  to  be  largely 
due  to  interpolation  ambiguities  and  not  a  fault  in  the  data. 

Contours  resulting  from  model  predictions  are  presented  in  Figures  4.3-4.6,4.9-4.12,  and  4.15- 
4.18.  These  contours  were  determined  from  predictions  for  a  rectangular  grid  (7x11)  and  do  not  appear 
to  be  in  especially  close  agreement  with  the  data.  The  general  features  and  order  of  magnitude  of 
predicted  and  observed  contours  appear  to  be  similar,  however.  Figures  4. 1,4.7,  and  4.13  present 
model/data  comparisons  at  the  actual  sampler  locations.  Those  figures  generally  reveal  a  fair 
representation  of  the  data.  Some  distortion  cannot  be  avoided  as  a  result  of  interpolation  ambiguities 
(especially  on  an  irregular  grid)  and  not  a  fault  of  the  data  themselves. 

4.2.  Discussion  of  the  Model/Data  Comparisons  for  the  Dugway  Tests 

4.  2  .1  Test  T0009 

For  this  test,  the  source  was  initially  located  at  J4  (0,0)  rather  than  (-4,30)  as  indicated  in  the  data 
figure.  However,  because  the  bulk  of  the  plume  missed  the  sampling  grid  to  the  west  (toward  the  bottom 
of  the  figure)  from  this  release  point,  the  source  was  relocated  to  (-4,30)  within  the  first  15  minutes  of  the 
test.  Because  the  smoke  release  continued  during  the  relocation,  the  data  closest  to  the  release  point 
appears  to  suggest  a  release  point  between  (0,0)  and  (-4,30). 

Given  this  complication,  one  may  make  the  following  conclusions  from  examining  Figures  4.1  to  4.6: 

(a)  Agreement  between  models  and  data  are  qualitatively  good, 

(b)  Model  predictions  are  generally  within  a  factor  of  3  of  the  data.  The  average  concentration  data 
have  larger  values  than  the  model  predictions  for  the  25  m  transect,  but  become  in  closer  agreement  with 
the  model  predictions  for  the  further  transects  50-200  m. 

(c)  As  can  best  be  seen  with  the  limited  data  available  (especially  at  the  further  transects),  the  lateral 
spreading  of  the  plume  seems  to  be  reasonably  well  predicted  by  the  models. 

(d)  In  terms  of  contour  shapes  and  magnitudes,  the  models  seem  to  be  in  reasonable  agreement 
with  the  data.  The  RIMPUFF  model  predicts  irregular  looking  contours  perhaps  because  of  the  very  short 
meteorological  time  step  used  in  that  model  and  thus  the  many  different  directions  the  time-dependent 
puff  releases  follow  during  the  experiment. 

4.  2  .2  Test  TO0 10 

The  wind  direction  frequency  was  found  to  be  bimodal  for  this  test  whereas  the  wind  direction 
frequency  distribution  for  tests  T0009  and  T0011  were  both  nearly  normal.  A  mode!  that  ’.ccounts  for 
time-dependent  variations  in  wind  direction  is  really  needed  in  this  case  because  a  single  average  wind 
direction  does  not  fully  represent  the  distribution  in  wind  directions  over  the  test  period. 

The  following  conclusions  may  be  drawn  from  the  examination  of  Figures  4.7  to  4.12. 
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(a)  Model  predictions  are  within  a  factor-of-2  (model  predictions  are  larger)  of  the  data  for  the  25  m 
and  50  m  transects.  However,  model  predictions  are  significantly  higher  than  the  data  for  the  100  m  and 
200  m  transects.  A  careful  examination  of  the  frequency  distribution  and  the  contour  plots  for  T0010 
shows  that  the  data  does  reflect  the  observed  wind  pattern.  The  plume  passing  off  the  sampling  grid 
makes  it  difficult  to  evaluate  the  cause  of  the  discrepancies  at  the  100  m  and  200  m  arc. 

(b)  Lateral  spreading  appears  to  be  consistent  between  the  models  and  data  --  given  the  limited  data 
available  on  the  100  m  and  200  m  transects. 

(c)  Concentration  maxima  arc  at  approximate!*  i'ie  correct  location  on  the  25  and  50  m  transects.  No 
evaluation  can  be  made  for  the  100  and  200m  transects  due  to  lack  of  a  refined  grid  of  data  values,  and 

(d)  The  predicted  average  concentration  contours  are  good  representations  of  the  data  contour 
with  the  RIMPUFF  model  providing  puffs  that  are  more  expansive  laterally  and  longitudinally  as  compared 
with  the  data. 

4.2.3  Test  TO0 11 

During  this  test,  the  wind  direction  changed  while  the  generators  were  being  field-repaired.  To 
accommodate  this  wind  shift,  the  generator  was  relocated  prior  to  restarting  the  test.  Both  locations  are 
indicated  in  Figure  4.14.  Conclusions  from  the  model/data  comparisons  for  test  T001 1  ere  as  follows: 

(a)  Model  predictions  are  excellent  for  this  case  both  in  terms  of  peak  concentration  value  at  each 
transect  and  lateral  spreading  at  each  transect.  Model  predictions  are  within  a  factor-of-2  of  the  data  at 
each  transect, 

(b)  Movement  of  the  plume  off  the  grid  at  the  200  m  transect  prevents  a  detailed  comparison  at  the 
distance, 

(c)  All  model  predicted  contours  provide  reasonable  representations  of  the  data  with  the  model 
predictions  showing  less  meandering  than  the  data.  However,  the  distortion  in  the  data  contour  may  be 
largely  due  to  the  attempt  to  interpolate  data  values  on  an  irregular  and  somewhat  sparse  grid. 

4.3  Concluding  Remarks  --  Dugway  Tests 

One  may  conclude  from  examination  of  the  model/data  comparisons  that,  for  all  nonzero  data  values 
and  all  the  models  combined,  the  model  predictions  give  50%  within  a  factor  of  two  and  70%  within  a 
factor  of  three  of  the  data  The  Dugway  field  studies  data  base  does  not  contain  a  large  enough  set  of  data 
for  a  definitive  evaluation  of  models  under  flat  terrain  conditions.  In  addition,  data  were  collected  only 
under  unstable  meteorological  conditions  and  short  distances  downwind  (up  to  200  m).  No  data  for 
neutral  (except  for  the  first  part  of  test  T001 1  )  or  stable  meteorological  conditions  were  obtained. 

The  performance  of  each  of  the  models  (under  these  unstable  conditions  and  short  downwind 
distances)  appeared  to  be  quite  reasonable,  however,  considering  the  simplicity  of  the  theories  and  the 
limited  data  required  for  model  input. 

The  presence  of  a  rising  centerline  observed  beyond  about  200  m  provides  a  unique  challenge  to 
the  models,  but  that  phenomenon  was  not  measured  quantitatively  in  terms  of  average  concentrations 
during  the  Dugway  tests.  The  rising  centerline  phenomenon  was  observed  and  quantified  as  part  of  the 
Camp  Atterbury  tests  and  the  performance  of  the  models  is  presented  later. 

A  final  comparison  of  the  Dugway  average  concentration  data  and  the  model  predictions  are  given  in 
Figure  4.19.  This  figure  is  actually  a  combined  plot  of  results  presented  in  Figures  4.1 , 4.7,  and  4.13.  A  45 
degree  line  on  the  plot  (not  shown)  would  indicate  perfect  agreement  between  predictions  and 
observations.  The  two  factor-of-2  lines  represent  predictions  that  are  a  factor  of  two  greater  or  less  than 
the  observed  value.  The  two  factor-of-3  'mes  on  the  plot  represent  predictions  that  are  a  factor  of  three 
greater  or  less  than  the  observed  value.  Note  the  large  number  of  points  within  the  factor-of-3  envelope. 
Points  on  one  axis  or  the  other  indicate  spatial  offset  between  the  observed  and  predicted  plumes  (i.e., 
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Figure  4.19  Scatterplot  comparison  of  predicted  and  observed  average 
concentrations  for  all  four  models  to  field  data  acquired  for  all 
three  Dugway  field  tests. 
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the  prediction  is  zero  and  the  associated  observation  is  nonzero  or  vice-versa).  The  typical  range  of 
precision  believed  representative  of  Gaussian  model  evaluations  is  a  factor  of  two  for  50%  of  the  time.  It  is 
now  said  for  the  first  time  that  the  particulate  phase  of  the  fog  oil  plume  (for  these  short  distances  and 
stability  classes  B,  C,  and  D)  can  be  represented  by  the  Gaussian  dispersion  model. 

5.0  THE  CAMP  ATTERBURY  FIELD  STUDIES 

This  section  provides  a  short  summary  of  the  data  acquired  in  the  October  /  November  1987  field  tests  at 
Camp  Atterbury  (Liljegren  et  al.,  1989).  Four  data  sets  are  used  for  model  testing  in  the  next  section.  More 
deiail  concerning  all  the  data  acquired  as  well  as  the  equipment  and  methodology  used  is  given  in  a  data 
report  by  Liljegren  (Liljegren  et  al.,  1989). 

5.1  The  Site  and  the  Sampling  Layout 

A  field  study  on  fog-oil  dispersion  was  carried  out  over  a  five-week  period  in  October  and  November 
1987  at  the  Atterbury  Reserve  Forces  Training  Center  (ARFTC)  near  Columbus,  Indiana.  The  test  site  was 
located  in  area  7A  in  the  ARFTC  complex.  The  topography  of  the  site  and  surrounding  terrain  are 
presented  in  Fig.  5.1.  The  site  is  positioned  in  a  meadow  among  low  hills  which  rise  approximately  50  m 
above  the  mean  site  elevation  of  218  m  above  mean  sea  level.  The  meadow  was  covered  primarily  with 
grass  about  1  m  high.  The  lower  portions  of  the  site,  where  water  collected  after  a  rainstorm,  were  covered 
with  briars,  brambles  and  reeds  1  to  2  m  high  in  addition  to  grass.  The  surrounding  hills  were  densely 
wooded  with  deciduous  trees  approximately  20  m  high.  Prior  to  the  testing  period,  which  spanned  the  first 
two  weeks  of  November,  the  trees  had  dropped  most  of  their  leaves. 

The  sampling  network  used  at  Camp  Atterbury  is  illustrated  schematically  in  Fig.  5.2.  The  filled 
circles  indicate  the  positions  of  the  masts  supporting  the  concentration  samplers.  The  network  was 
composed  of  50  sampling  locations  on  five  linear  transects  oriented  perpendicular  to  the  anticipated 
prevailing  wind  direction  (225  deg  east  of  north).  In  order  to  resolve  any  near-ground  concentration 
gradients,  concentration  samplers  were  deployed  at  elevations  of  1 ,2,4  and  8  m  on  the  first  four  transects 
and  at  elevations  of  2  and  8  m  on  Transect  5  for  a  total  of  192  samplers  per  test.  The  transects  were 
arranged  at  distances  of  50,  100,  250,450  and  675  m  from  the  baseline.  The  fifth  transect  was  located  10 
to  25  m  upwind  of  a  dense  line  of  trees  separating  Area  7A  and  Area  1A.  Because  the  plume  was 
expected  to  spread  linearly  as  it  traveled  downwind  from  the  source,  the  sampler  spacing  increased  on 
each  transect  in  order  to  maintain  approximately  constant  lateral  resolution.  The  samplers  were  spaced  at 
intervals  of  15.2,30.4,45.7,61.0,  and  122  m  respectively  on  transects  1  to  5. 

A  number  of  possible  source  locations  were  laid  out  depending  upon  wind  direction.  Real-time 
plume  model  predictions  made  with  a  computer  in  the  field  were  used  to  locate  the  source  for  each  test 
depending  upon  wind  direction  and  the  orientation  of  sampling  transects.  In  this  way,  efforts  were  made  to 
have  the  plume  from  each  test  cover  as  much  of  the  sampling  grid  as  possible.  The  location  of  the  source 
and  prevailing  wind  direction  for  the  four  fog-oil  tests  are  presented  in  Fig.  5.3a  and  5.3b.  The  tests  are 
designated  by  the  date  and  test  number  on  that  date.  For  example,  the  third  fog  oil  test,  designated 
1104872,  was  the  second  test  executed  on  November  4,  1987. 

No  on-site  meteorological  data  were  available,  and  the  closest  National  Weather  Service  Station  is  in 
Indianapolis,  Indiana.  Other  than  assisting  in  the  determination  of  the  possible  best  times  during  the  day 
for  a  field  test,  such  off-site  data  was  of  little  use  considering  that  on-site  meteorological  data  from  a  10  m 
tower  was  acquired  at  the  site  for  each  of  the  experiments. 

5.2.  Data  on  the  Fog-Oil  Smoke  Source 

A  M3A4  military  smoke  generator  was  used  to  produce  the  smoke  by  vaporizing  fog  oil  and  ejecting 
it  into  the  atmosphere  at  a  nominal  rate  of  46  g/s  according  to  the  M3A4  Operator's  Manual.  The  resulting 
aerosol,  as  with  the  M3A3E3  smoke  generator,  is  a  two-phase  mixture  of  droplets  and  vapor.  The 
operation  of  a  single  M3A4  smoke  generator  from  a  release  height  of  2  m  provides  a  close  approximation 
to  a  continuous,  ground-level  source  of  varying  emission  rate  with  time.  Averages  of  the  fog  oil  source 
data  are  presented  in  Table  5.1.  During  test  1103871  the  power  supply  for  the  portable  computer  used  to 
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Figure  5.1  Topological  Map  of  the  Camp  Atterbury  Dispersion  Site.  Elevations  are  in 
feet  above  sea  level  with  isopleths  in  increments  of  10  feet.  The  horizontal 
grid  is  in  Universal  Transverse  Mercator  coordinates,  with  marked 
increments  of  1  kilometer.  Sampling  mast  locations  are  indicated  with  filled 
squares!  a  ].  Circles  [  °|  indicate  the  various  release  points.  Particle  size 
measurements  were  carried  out  at  the  locations  marked  by  diamondsfo].  The 
location  of  the  10  m  meteorological  instrument  tower  is  indicated  by  the 
filled  circlef  •].  The  locations  of  the  2  m  wind  monitoring  stations  operated 
prior  to  the  dispersion  tests  for  the  purpose  of  assessing  the 
micrometeorology  of  the  site  are  indicated  by  the  filled  diamondl  ♦]• 
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Figure  5.2  Nominal  sampling  network  for  Atterbury  87  dispersion  field  study 
Symbols  are  the  same  as  for  Fig.  5.1. 
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Table  5.1  Summary  of  source  data  for  fog  oil  tests  conducted  during  Atterbury-8/. 


test: 

1 1 03871 T 

1104871 

1104872 

1  106871 

Source  Location  (m) 

CO 

b 

l 

3.0,  -4.7 

-47.7,  47.7 

-47.7,  47.7 

Mass  Released  (kg) 

115.9 

44.0 

1 17.8 

193.0 

Release  Duration  (min) 

55.9 

27.5 

48.2 

76.2 

Rolenne  Rail.'  (g/s) 

34  .0 

26.7 

4  0.7 

42.2 

Exit  Temperature  (°C) 

486 

444 

401 

415 

Wind  Speed  (m/s) 

2.6 

2.3 

1.7 

0.8 

Wind  Direction  (°E  of  N) 

241 

249 

247 

224 

09  n 

28.6 

23.4 

26.0 

33.8 

Temperature  (°C) 

23.4 

20.0 

25.3 

8.0 

Relative  Humidity  (%) 

41 

25 

43 

4  0 

^Based  on  last  live  minutes  of  source  data  record  due  to  equipment  malfunction. 
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Figure  5.4  One-minute  averaged  source  data  for  test  1104871  at  Camp  Atterbury. 

Top:wind  speed  and  wind  direction;  bottom:exit  port  temperatures  and 
release  rate. 


nd  Direction  (°E  of  N) 


Temperature  (°C) 


Figure  5.5A 


\j  4'i  ~A//\  aaa 

V  « 


TEST:  1104872 

SMOKE:  FOG  OIL 
START:  15:26 
END:  16:14 

-  Speed 

-  - - Direction 


i  -  /. 

s  1  v  > 


Figure  5.5B 


Time  Relative  to  Start  (minutes) 


TEST: 

SMOKE; 

START: 

END: 


1104872 
FOG  OIL 
15:26 
16:14 

Tcvrr  . 


'Exrr.3 

Rate 


,r\\ 


Time  Relative  to  Start  (minutes) 

Figure  5.5  One-minute  averaged  source  data  for  test  1104872  at  Camp  Atterbury. 

Top:wind  speed  and  wind  direction;  bottom:exit  port  temperatures  and 
release  rate. 
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Figure  5.6 


One-minute  averaged  source  data  for  test  1 106871  at  Camp  Atterbury. 
Top:wind  speed  and  wind  direction;  bottonrexit  port  temperatures  and 
release  rate. 


record  the  digital  output  of  the  source  data  logger  failed.  As  a  result  only  that  portion  of  the  data  record 
covering  the  last  five  minutes  of  the  test  could  be  recovered  from  the  datalogger.  Because  the  test  was 
ultimately  terminated  by  a  failure  of  the  M3A4  smoke  generator  unit,  the  average  exit  temperature  during 
the  last  five  minutes  of  the  test  may  not  be  representative  of  the  entire  test. 

Plots  of  the  1 -minute  averages  of  wind  speed,  wind  direction,  exit  temperature  and  release  rate  are 
presented  in  Figs  5.4  to  5.6  for  tests  1104871,  1104872  and  1106871  respectively.  For  test  1104871, 
Fig.  5.4a  reveals  that  the  wind  speed  was  highly  variable  although  the  wind  direction  was  fairly  constant. 
As  a  result,  the  material  plume  tended  to  move  downwind  in  pulses  rather  than  as  a  uniform  cloud.  Fig. 
5.4b  confirms  the  previous  statement  that  the  M3A4  failed  repeatedly  during  this  test.  Once  the  pulse-jet 
engine  was  started,  the  exit  temperature  climbed  rapidly  until  the  flow  of  fog  oil  was  initiated.  As  the  oil  flow 
increased,  the  release  rate  increased  as  well  and  the  exit  temperature  declined  to  a  steady  state  value  of 
approximately  450  deg  C.  About  18  minutes  into  the  test,  the  pulse-jet  engine  failed  for  the  first  time,  as 
indicated  by  the  precipitous  drop  in  exit  temperature  and  release  rate.  The  pulse-jet  engine  was  restarted 
within  one  minute  and  the  test  continued.  This  pattern  was  repeated  until  after  the  third  failure  of  the 
pulse-jet  engine,  whereupon  the  test  was  terminated. 

As  Fig.  5.5a  and  5.5b  indicates,  test  1104872  was  more  successful.  The  wind  speed  and  direction 
were  very  similar  to  the  earlier  test  but  the  M3A4  performed  much  better.  As  before,  the  exit  temperature 
shows  an  initial  "spike"  as  the  pulse-jet  engine  was  started  and  allowed  to  warm  up.  Oil  flow  was  initiated 
more  slowly  but  increased  gradually  throughout  the  test,  as  indicated  by  the  gradual  decline  in  the  exit 
temperatures.  The  variation  in  the  measured  exit  temperatures  among  the  three  exhaust  ports  is  due  to 
small  variations  in  the  proximity  of  the  thermocouples'  measuring  junction  to  the  exit  plane.  The  large 
dynamic  pressure  on  the  thermocouples  due  to  the  high  exit  velocity  caused  them  to  bend  away  from  the 
exit  plane.  As  the  exhaust  temperature  declined  rapidly  with  distance  from  the  ports,  small  variations  in  the 
thermocouple  positions  produced  significant  variations  in  the  measured  temperatures. 

The  smoke  generator  performed  flawlessly  during  test  1106871  (see  Figures  5.6a  and  5.6b).  As 
Fig.  5.6b  indicates,  the  wind  speed  was  much  lower  than  for  the  previous  tests  and  the  wind  direction  was 
highly  variable  and  resulted  in  a  much  larger  value  for  (Je  in  Table  5.1.  Despite  the  variability  of  the  wind 
direction,  the  smoke  plume  remained  on  the  sampling  network.  Fig.  5.6b  reveals  that  oil  flow  was  initiated 
very  soon  after  starting  the  pulse-jet  engine  and  the  exit  temperature  stabilized  rapidly.  The  data  also 
reveal  that  release  rate  was  much  more  variable  during  the  test.  This  is  due  to  the  lower  ambient 
temperature  (8  deg  C  vs.  25  deg.  C)  which  correspondingly  increased  the  viscosity  of  the  fog  oil.  This  in 
turn  affected  the  ability  of  the  fog  oil  pump  to  draw  the  oil  from  the  supply  drum  and  force  it  through  a  spray 
nozzle  into  the  smoke  generator.  Interestingly,  this  does  not  seem  to  have  had  much  effect  on  either  the 
overall  release  rate  or  the  overall  exit  temperature  presented  in  Table  5.1 . 

Droplet-size  data  acquired  during  these  fog  oil  tests  show  that  the  geometric  mass  mean  diameter  of 
the  fog  oil  aerosol  is  about  1  micron.  This  implies  that  no  deposition  is  expected  (also  verified  in  the 
Dugway  and  Camp  Atterbury  field  tests);  further,  laboratory  and  field  studies  verify  that  no  significant 
evaporation  will  occur  over  the  distances  that  concentration  data  were  obtained.  Finally,  the  Dugway  field 
tests  and  supporting  laboratory  studies  done  at  the  University  of  Illinois  indicate  that  there  is  no  significant 
vapor  phase  (less  than  1%  of  mass  emitted)  to  the  fog  oil  emission.  Consequently,  the  fog  oil  should  act 
like  a  tracer  with  99%  contained  in  the  aerosol  phase. 

6.0  COMPARISON  OF  MODELS  WITH  FIELD  DATA  TAKEN  AT  CAMP 

ATTERBURY 


6.1  Preparation  of  Model  Inputs 

As  in  the  Dugway  field  tests,  the  data  from  the  source  and  meteorological  towers  were  used  to 
determine  the  required  inputs  to  the  four  models.  Table  6.1  preserts  a  summary  of  the  main  source  and 
meteorological  variables  for  each  of  the  four  Camp  Atterbury  tests.  As  in  the  Dugway  tests,  all  the  models 
used  10  meter  level  wind  speeds  and  directions,  which  were  chosen  based  on  discussions  with  the 
modelers.  Pasquill  stability  classes  are  used  in  two  models:  BEAR  and  INPUFF-PG.  Stability  classes  were 
all  B  except  for  case  1106871  which  had  stability  class  A.  The  RIMPUFF  and  INPUFF-Onsite  models  used 
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Test  Designation 

1103871 

1104871 

1104872 

1106871 

Smoke  Type 

Fog  Oil 

Fog  Oil 

Fog  Oil 

Fog  Oil 

Begin  Test  (CDT) 

10:31:06 

09:36:33 

15:25:50 

10:51:50 

End  Test  (CDT) 

11:27:00 

10:06:00 

16:14:00 

12:08:00 

Duration 

00:55 

00:29 

00:48 

01:16 

Transects  in  Operation 

1|2|3|4|5 

1|2|3|4|5 

1|2|3|4|5 

1|2|3|4|5 

Source 

Coordinates{x(m)  ,y(m)} 

3.0, -4.7 

3.0, -4.7 

-47.4,47.4 

-47.4,47.4 

Mass  Released  (kg) 

115.9 

44.0 

117.8 

193.0 

Release  Time  (min) 

55.9 

27.5 

48.2 

76.2 

Rate  of  Release  (g/s) 

34.6 

26.7 

40.7 

42.2 

2  m  Wind  Velocity  (m/s) 

2.6 

2.3 

1.7 

0.8 

2  m  Wind  Direction  (°E  of  N) 

241 

249 

247 

224 

2  m  Sigma-theta  (°) 

28.6 

23.4 

26.0 

33.8 

2  m  Ambient  Temp  (°C) 

23.4 

20.0 

25.3 

8.0 

2  m  Relative  Humidity  (%) 

41 

25 

43 

40 

Particle  Size _ 

Location[x(m),y(m)] 

316.4,298.4 

316.4,298.4 

316.4,298.4 

184.0,164.5 

Mass  Mean  Dia.  (pm) 

0.96 

- 

0.67 

1.01 

Geometric  Mean  Dia.  (pm) 

0.86 

- 

0.56 

0.90 

Geometric  Std.  Deviation 

1.55 

• 

1.71 

1.53 

Meteorology _ 

10  m  Wind  Velocity  (m/s) 

5.5 

5.1 

4.7 

1.6 

1 0  m  Wind  Direction  (°E  of  N) 

239 

249 

261 

240 

10  m  Wind  Inclination  (°) 

-10 

-12 

-11 

-8 

10  m  Sigma-theta  (°) 

16.2 

11.0 

18.6 

35.4 

10  m  Sigma-phi  (°) 

8.8 

8.2 

8.2 

14.9 

10  m  Sigma-u  (m/s) 

1.52 

1.17 

1.26 

0.70 

10  m  Sigma-v  (m/s) 

1.49 

0.94 

1.54 

0.90 

10  m  Sigma-w  (m/s) 

0.94 

0.82 

0.78 

0.49 

1 0  m  •  2  m  AT  (°C) 

-0.41 

- 

1.89 

-0.89 

Wind  Power  Law  Exponent 

0.143 

0.149 

0.183 

0.103 

Est.  Obukhov  Lenght  (m) 

-5 

-6 

-3 

-5 

Est.  Friction  Velocity  (m/s) 

0.909 

0.819 

0.871 

0.271 

Est.  Roughness  Height  (m) 

0.20 

0.20 

0.20 

0.20 

Est.  Inversion  Height  (m) 

700 

700 

700 

700 

Est.  Convection  Velocity 

2.53 

2.25 

2.07 

1.29 

zi/L  =  -k(w7u'  )A3 

-9.93 

-8.53 

-5.52 

-44.00 

Notes _ 

1 .  All  values  are  averaged  over  the  period  BEGIN-10  to  END  +10  minutes. 

2.  Only  last  9  minutes  of  source  data  record  is  available  for  tests  1103871  &  1113871. 

3.  Only  20  minutes  of  met  record  is  available  for  test  1104871. 

4.  No  particle  size  measurements  were  obtained  during  test  1 1 04871  due  to  insufficient  smoke 
concentrations  at  the  instrument  location. 

5.  Convection  velocity  from  Deardorff  and  Willis  (1975)  correlation 

6.  Negative  mean  wind  inclinations  reflect  the  mean  subsidence  flow  observed  under  convective 
conditions  (Liljegren  1989) 


Table  6.1  Summary  of  data  for  four  field  tests  at  Camp  Atterbury  for  fog-oil  dispersion. 
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Comparison  of  RIMPUFF,  BEAR,  INPUFF-Onsite  and  INPUFF-PG 
mathematical  models  with  average  concentrations  from  test  1103871  at 
Camp  Atterbury. 


Figure  6.3.  Comparison  of  RIMPUFF,  BEAR,  INPUFF-Onsite  and  INPUFF-PG 

mathematical  models  with  average  concentrations  from  test  1 104872  at 
Camp  Atterbury. 


6.4.  Comparison  of  RIMPUFF,  BEAR,  INPUFF-Onsite  and  INPUFF  PG 

malhemalical  models  with  average  concentrations  (rom  test  1 106871  at 
Camp  Atterbury. 


Ce  and  a<s  values  from  the  10-m  level  as  a  function  of  time.  As  with  the  application  of  the  models  to  the 
Dugway  data,  the  RIMPUFF  model  used  15  second  averaged  meteorological  data  and  the  other  three 
models  used  3  minute  averaged  meteorological  data.  That  data  was  available  and  provided  to  the  model. 

The  release  height  was  chosen  to  be  2  m  above  the  ground  for  all  models  since  it  represented  the 
observed  height  above  ground  of  the  center  of  the  initial  puff.  The  strong  initial  mixing  that  takes  place 
and  the  slight  buoyancy  of  the  plume  at  exit  leads  to  some  uncertainty  as  to  the  effective  height  of 
release. 

The  release  rate  of  fog  oil  as  a  function  of  time  from  the  smoke  generator  was  measured  and  was 
input  to  each  of  the  four  models.  The  exit  temperature  as  a  function  of  time  was  not  used  in  the  models 
since  they  could  not  use  that  information.  The  assumption  of  a  nonbuoyant  source,  in  the  application  of 
these  models,  precludes  the  use  of  those  exit  temperature  data.  In  essence,  the  models  are  not  complex 
enough  in  structure  to  require  some  of  the  measured  data  as  input  to  the  codes.  More  complex  models 
may  be  tested  in  the  future  which  would  use  all  of  the  data  measured.  In  any  case,  the  data  collected  were 
necessary  to  develop  a  good  understanding  of  the  physics  of  the  fog  oil  plume  dispersion. 

6.2  Results  of  the  Model/Data  Comparisons 

Figures  6. 1-6.4  present  the  comparisons  of  average  concentrations  with  lateral  distance  for  each  of 
the  four  models  at  each  of  the  five  downwind  transects.  The  following  discussion  can  be  made  of  those 
comparisons: 

(a)  Model  agreement  with  the  data  is  much  better  close  to  the  source  (at  250  m  from  the  generator 
and  closer)  than  further  from  the  source  (at  the  450  and  675  m  transects).  The  rise  of  the  plume  centerline 
noticed  beyond  about  250  m  leads  to  much  smaller  concentrations  measured  than  predicted.  Model 
predictions  assume  that  all  puffs  released  will  remain  at  the  height  of  release  at  all  distances  downwind  with 
no  convective  rise.  The  general  overprediction  of  the  models  at  the  450  m  and  675  m  transects  is  likely 
due  to  the  lack  of  inclusion  of  convective  rise  of  the  puffs  released  from  the  generator, 

(b)  The  BEAR  and  INPUFF-PG  models  are  nearly  identical  in  predictions  for  all  cases.  The  theory  of 
the  models  is  very  similar  (both  use  the  Pasquill-Giflord  stability  class  method)  and  both  use  time- 
dependent  winds  at  the  10-m  level  to  transport  those  dispersing  puffs, 

(c)  The  RIMPUFF  model  predicts  higher  average  concentrations  than  the  other  models  for  most 
longitunal  distances, 

(d)  The  BEAR,  INPUFF-Onsite,  and  INPUFF-PG  models  predict  closer  to  each  other  than  to  the 
data.  One  may  conclude  from  examination  of  the  model/data  comparisons  that,  for  nonzero  data  values 
and  all  the  models  combined,  the  model  predictions  give  32%  within  a  factor  of  two  and  46%  within  a 
factor  of  three  of  the  data  . 

One  troubling  aspect  to  the  model/data  comparisons  is  the  generally  poor  agreement  in  the  rate  of 
concentration  decay  with  distance  between  model  predictions  and  data.  At  distances  between  about  50 
and  250  m  the  model  predictions  are  generally  within  a  factor  of  two  to  three.  However,  for  distances  of  50 
m  or  less,  the  models  tend  to  significantly  underpredict  or,  in  some  cases,  overpredict  the  data.  At 
distances  beyond  250  m,  the  models  tend  to  significantly  overpredict  the  data  and  this  disagreement  is 
due,  in  part,  to  the  rising  centerline  phenomenon. 

The  Gaussian  puff  models  suffer  from  the  assumption  that  the  puff  centerlines  must  always  stay  at 
the  height  of  release  and  do  not  rise.  There  appears  to  be  no  judicious  choice  of  the  oz  parameter  of  the 
Gaussian  model  to  "correct"  for  the  error  of  not  allowing  for  a  rising  centerline.  The  data  showed  that  the 
lateral  distributions  of  concentration  were  approximately  Gaussian  in  shape.  The  poor  agreement  between 
Gaussian  models  and  the  data  for  the  rate  of  centerline  concentration  decay  with  distance  underscores 
the  fact  that  the  vertical  profile  of  concentration  is  not  Gaussian  in  shape. 

A  final  comparison  of  the  Camp  Atterbury  average  concentration  data  and  the  model  predictions  are 
given  in  Figure  6.5.  This  figure  is  actually  a  combined  plot  of  results  presented  in  Figures  6.1  to  6.4.  The 
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Figure  6.5  Scatterplot  comparison  of  predicted  and  observed  average 
concentrations  for  all  four  models  to  field  data  acquired  for  all 
four  Camp  Atterbury  field  tests. 
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45-degree  line  on  the  plot  (not  shown)  represents  perfect  agreement  between  predictions  and 
observations.  The  two  factor-of-2  and  the  two  factor-of-3  lines  represent  predictions  that  are  a  factor  of 
two  and  factor  of  three  greater  or  less  than  the  observed  values,  respectively.  Note,  as  in  the  Dugway 
tests,  the  large  number  of  points  within  the  factor-of-3  envelope.  As  in  the  Dugway  comparison,  many 
predicted/observed  points  are  close  to  the  origin.  The  typical  range  of  precision  believed  representative 
of  Gaussian  model  evaluations  is  a  factor  of  two  for  50%  of  the  time.  It  can  be  said  at  this  site,  just  like 
Dugway,  that  the  particulate  phase  of  the  fog  oil  plume  (for  these  short  distances  and  stability  classes  B,  C, 
and  D)  can  be  represented  by  the  Gaussian  dispersion  model. 

7.0  SUMMARY  AND  CONCLUSIONS 

Seven  data  sets  for  fog-oil  smoke  dispersion  in  flat  terrain  have  been  used  to  test  four  Gaussian  puff 
models.  The  models  are  simple  in  formulation  and  require  time-dependent  source  and  meteorological 
data  for  input.  The  models  divide  the  release  into  a  series  of  puffs  which  are  dispersed  and  transported  in 
the  local  wind  direction.  The  puff  centerlines  remain  at  the  height  of  release  during  their  lifetimes  but  the 
puffs  travel  in  directions  that  change  with  the  wind  direction  data.  The  models  differ  in  terms  of  their 
meteorological  time-averaging  period  and  their  treatment  of  puff  dispersion:  use  of  Pasquill-Gifford 
stability  classes,  onsite  turbulence  data  to  provide  direct  estimates  of  dispersion  coefficients  (avoiding  the 
P-G  stability  class  scheme),  or  use  of  the  Smith-Hay  formulas  for  puff  dispersion  (also  avoiding  the  P-G 
stability  class  scheme). 

The  data  sets  represent  3  data  cases  from  the  experiments  at  Dugway  Proving  Ground  and  4  data 
cases  from  the  Camp  Atterbury  experiments.  Data  were  measured  at  several  transects  downwind  in  each 
experiment  with  the  Dugway  data  obtained  25  m  -  200  m  downwind  and  the  Camp  Atterbury  data  obtained 
from  50  m  -  675  m  downwind.  Smoke  generator  runs  varied  from  29-98  minutes.  The  concentration  data 
represented  average  values  over  the  period  of  the  experiment.  In  the  Camp  Atterbury  experiments,  data 
were  obtained  at  four  levels  from  1  m  to  8  m  .  The  Dugway  concentration  data  were  obtained  a*  ?  m  above 
the  ground  at  all  sampling  stations. 

The  results  show  that  the  four  Gaussian  puff  models  tested  can  predict  within  a  factor  of  three  under 
the  convective  and  neutral  conditions  tested  to  distances  of  about  250  m.  Beyond  that  distance,  the 
plume  tends  to  rise  leading  the  models  to  significantly  overpredict  average  concentrations  at  ground  level. 
No  data  were  acquired  at  these  two  sites  under  stable  conditions.  (Fog-oil  plumes  under  stable  conditions 
were  measured  at  the  Meadowbrook  site  in  Red  Bluff,  California  but  model  testing  with  those  data  is  not 
the  subject  of  this  report. ) 

The  factor  of  two  to  three  accuracy  for  this  short-distance  dispersion  has  been  seen  consistently  for 
over  the  past  twenty  years  with  respect  to  point  source  tracer  releases  near  the  ground.  The  failure  of  the 
models  to  predict  plume  rise  from  the  ground  under  convective  conditions  beyond  about  250  m  is  a  major 
flaw  which  needs  to  be  emphasized.  Given  the  accuracy  of  field  data  generally,  the  evaluation  of  the 
Gaussian  models  will  yield  similar  results  consistently,  barring  the  possibility  of  poor  source 
characterization.  The  present  work  adds  to  the  state-of-the-art  in  that  it  is  now  shown  that  the  particulate 
phase  of  a  fog  oil  plume  acts  like  a  tracer  in  its  dispersion  in  the  atmosphere  -  for  the  short  distances  and 
stability  classes  tested  with  the  Dugway  and  Camp  Atterbury  data. 

Although  the  Gaussian  models  tested  with  the  fog-oil  data  indicated  general  agreement  with  the 
average  concentration  data  within  a  factor  of  two  or  three  for  the  average  concentrations  measured  (for 
about  50%  of  data  points),  there  were  problems  with  the  models  (under  these  largely  unstable 
atmospheric  conditions)  in  that 

(i)  such  agreement  was  only  within  the  first  250  m  from  the  smoke  generator  with  more  significant 
discrepancies  at  further  distances, 

(ii)  there  is  a  tendency  to  significantly  overpredict  or  underpredict  average  concentrations  at 
distances  less  than  100  m  and  to  significantly  overpredict  at  distances  greater  than  about  250  m.  The 
discrepancies  found  at  longer  distances  (greater  than  about  250  m)  is  due  to  a  rise  in  the  plume 
centerline,  and 
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(iii)  the  decay  of  concentration  with  distance  as  predicted  by  the  models  does  not  agree  with  the 
data  and  may  be  assumed  to  be,  in  part,  the  result  of  the  rising  centerline  seen  in  the  observed  plume. 
The  rising  centerline  phenomenon  has  been  observed  in  other  field  and  laboratory  studies  and  cannot  be 
represented  by  these  Gaussian  models.  A  new  stochastic  model  tested  with  the  data  seems  to  replicate 
the  systematic  behavior  in  the  data  (including  the  rising  centerline)  quite  well  (Liljegren,  1989). 

In  a  companion  report  by  Liljegren  (1989),  a  stochastic  model  is  compared  to  the  Camp  Atterbury 
data  and  is  shown  to  predict  the  rise  of  the  plume  under  convective  conditions  quite  well.  The  lack  of 
treatment  of  convective  turbulence  in  the  Gaussian  puff  models  is  apparently  the  cause  of  the 
discrepancies  between  the  model  predictions  and  the  data  beyond  about  250  m.  Future  work  would 
involve  the  testing  of  complex  terrain  versions  of  all  these  models  with  both  the  stable  and  unstable  plume 
dispersion  data  at  the  Meadowbrook  site. 

One  difficulty  with  the  application  of  the  four  Gaussian  models  to  the  seven  data  sets  is  the  simple 
nature  of  the  models  tested.  Indeed,  simple  models  are  needed  for  the  practical  and  rather  routine  smoke 
plume  modeling  used  by  the  U.S.  Army.  However,  simple  models  are  needed  for  some  of  the  practical  and 
rather  routine  smoke  plume  modeling  applications  used  by  The  U.S.  Army.  Some  of  the  questions  that  are 
raised  in  the  application  of  these  models  to  the  Dugway  and  Camp  Atterbury  data  sets  are: 

(a)  Should  the  wind  speed  at  the  release  height  be  used  or  is  the  10-m  value  more  representative? 
Or  is  an  average  over  the  height  of  the  plume  (at  a  given  transect)  more  representative9  An  entire  profile  of 
such  meteorological  data  is  available  from  the  experiments. 

(b)  What  is  the  initial  puff  radius  and  height  above  the  ground  once  the  initial  momentum  of  the 
plume  exiting  the  three  release  ports  is  accounted  for?  What  is  the  location  of  the  passive  diffusion  source 
(perhaps  located  upwind  of  the  smoke  generator)  that  is  "equivalent"  to  the  source  with  initial  momentum? 

(c)  What  is  the  stability  class  for  each  experiment?  The  uncertainty  in  estimating  the  stability  class  is 
obviously  present  and  the  results  are  definitely  dependent  or.  the  value  chosen. 

(d)  How  sensitive  are  the  predictions  to  the  kinds  of  uncertainties  that  have  been  identified  in  (a)  and 
(c)9 


There  are  no  definitive  answers  to  such  questions  (known  at  this  time)  and  the  models  are  clearly 
sensitive  to  the  choices  made.  In  any  case,  the  most  reasonable  choice  for  each  parameter  was  thought  to 
have  been  made  in  applying  the  models  to  the  data  in  this  report.  Furthermore,  the  agreement  between 
the  models  and  data  remains  reasonable  for  such  flat  terrain  applications  and  short  distances  for  which  the 
models  tested  well  (less  than  about  250  m).  The  accuracy  of  the  Gaussian  models  is  in  doubt  beyond  that 
distance  under  unstable  conditions.  The  problem  in  choosing  the  appropriate  parameters  for  model  input 
is  inherent  with  those  types  of  simple  Gaussian  puff  models. 
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B-02 


I 

During  these  Smoke  Weeks,  one  or  two  rows  of  aerosol  photometer  (AP)/che«1cal  Implnger  (Cl)  pairs 
H  were  placed  downwind  of  the  smoke  generator.  These  rows  were  approximately  normal  to  the  plume 

centerline  within  about  100  m  of  the  source.  Ho  data  exist  beyond  this  distance.  Aerosol  photometers 
I  measure  concentration  versus  time  through  an  optical  scattering  method.  These  instruments  are  fairly 

rapid  In  responding  and,  thus,  are  suitable  for  time-dependent  concentration  measurements,  but  their 

■  accuracy  In  terms  of  absolute  concentration  is  poor.  Thus,  the  AP  readings  were  calibrated  against 

measurements  of  dosage  made  simultaneously  by  the  chemical  Implnger.  i.e..  the  readings  were  scaled  on  a 

■  time-by-time  basis  such  that  concentration  dosages  are  equal  for  each  AP/CI  pair.  A  brief  description 
of  the  two  Smoke  Week  sites  and  available  data  follow. 

»  . . . . . .  _  „ . . .  ...  ... 

■|  sandy  soil,  scattered  pine  trees  and  brush  extended  to  an  elevation  of  16  m  above  ground  level.  Terrain 

9  In  the  vicinity  of  the  32-m  meteorological  tower  was  level  and  covered  with  grass,  except  for  rutted 

_  roadways.  Beyond  tne  roadways  there  was  a  margin  of  grass  with  brush  and  small  trees  ranging  from  0.5 

■  to  3.0  m  tall.  At  a  distance  of  300  m  to  the  north,  east,  and  south  there  was  a  thick  pine  forest 

_  15-20  m  tall  on  slightly  rolling  sand  hills.  The  roughness  length  was  estimated  to  be  50  cm  for  tne 

|  lower  levels  of  the  32-m  meteorological  tower  and  30  cm  for  the  upper  levels.  Roughness  determination 

was  complicated  by  wake  effects  of  obstacles  upwind. 

I 

Figures  2-4  present  sketches  of  the  layouts  of  the  smoke  generator,  meteorological  tower,  and  the 

■  aerosol  photometer/chemical  Implnger  system  for  the  Smoke  Wee*.  III  and  IV  trials.  The  main  sampling 
line  (LOS  2)  consisted  of  101  chemical  Implngers  (3  m  apart)  and  33  aerosol  photometers  (9  m  apart). 

M  The  auxiliary  sampling  line  (LOS  3)  consisted  of  27  chemical  Implngers  (6  m  apart)  and  9  aerosol 

photometers  (18  m  apart).  Model /data  comparisons  are  made  here  with  the  results  of  the  chemical 

■  Implngers.  Future  work  will  Involve  comparisons  of  models  on  a  concentration  versus  time  basis. 

■  Not  all  the  meteorological  data  measured  for  each  trial  could  be  used  by  the  models.  Available  was 

■ 

“  a  32-m  meteorological  tower  located  at  the  center  of  the  test  grid  which  was  equipped  to  measure 

™  (a)  horizontal  and  vertical  wind  directions  at  the  2.  4.  8.  16.  and  32-m  levels. 

m  (b)  horizontal  wind  speed  at  the  2.  4,  8.  16.  and  32-m  levels. 

1 

■  (c)  temperature  and  dew  point  at  the  2  and  32-m  levels,  and 

.  (d)  temperature  difference  between  the  0.5  and  4-m  levels. 
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ABSTRACT 


Numerous  mathematical  models  exist  to  predict  the  potential  environmental  Impact  resulting  from  the 
use  of  military  smokes.  Hazard  prediction  methods  require  the  modeling  of  smoke  dispersion  at  receptor 
locations  typically  20  to  2000  a  from  the  source.  Accurate  knowledge  of  smoke  concentrations  Is  also 
luportant  for  defining  toxicological  studies. 

This  paper  provides  a  preliminary  test  of  the  performance  of  four  mathematical  models  with  dosage 
data  from  Smoke  Week  III  Trials  4  and  19,  and  Smoke  Week  IV  Trial  3.  In  these  trials,  fog  oil  was 
released  from  a  H3A3  smoke  generator.  Although  these  data  were  not  obtained  for  the  explicit  purpose  of 
model  evaluation,  they  represent  the  best  data  available  at  the  present  time  for  this  purpose.  More 
extensive  data  are  being  obtained  as  part  of  the  present  project. 

The  Smoke  Week  data  used  for  model  testing  represent  fog-oil  trials.  In  a  typical  trial,  10  to  100 
dosage  samplers  are  set  out  along  one  or  two  transects  approximately  normal  to  the  plume  centerline  at 
distances  of  itcut  50-90  rs  downwind.  Meteorological  data  from  a  32-«  (45-o  In  Smoke  Week  IV)  Instrument 
tower  were  also  available. 

The  models  tested  are  C0MBIC,  ACT  II,  MAO  Puff,  and  the  Ludwig  (1977)  model.  All  models  are  of 
Gaussian  type.  C0MBIC.  ACT  II.  and  MAO  Puff  are  applied  to  meteorological  data  that  have  been  time- 
averaged  over  the  period  of  smoke  release.  The  Ludwig  Model,  on  the  other  hand.  Is  applied  to  time- 
varying  meteorological  data  over  the  trial  period.  C0MBIC  and  ACT  II  model  predictions  have  straight 
line  trajectories.  However,  the  MAD  Puff  "and  Ludwig  Models  have  variable  trajectories  due  to  the 
modeling  of  turbulence  along  the  trajectory  (for  MAO  Puff)  and  the  treatment  of  time-varying  wind 
direction  for  the  Ludwig  Model. 

Model/data  comparisons  reveal  that  all  four  models  are  within  factors  of  2-3  In  predicting  dosages 
at  each  of  the  rovs  of  samplers.  In  general,  the  models  tend  to  underestimate  lateral  spreading  and 
have  their  predicted  peak  values  offset  from  the  data.  It  was  found  that  for  such  short  release  times 
and  for  the  prediction  of  dosages,  the  lateral  movement  of  puff  centers  due  to  turbulence  (In  MAD  Puff) 
or  by  time-dependent  wind  direction  variations  (In  Ludwig)  do  not  appear  to  Improve  model  performance. 
Unfortunately,  the  rows  of  samplers  were  all  very  close  to  the  generator;  no  data  exist  beyond  about 
150  m  from  the  source.  Future  measurements  should  include  the  fog-oil  release  rate  from  the 
generator.  For  Smoke  Week  III,  Trials  4  and  19,  the  nominal  40  gal/hr  rate  of  release  had  to  be  assumed 
since  no  measurements  of  release  rate  or  total  fog  oil  consumed  were  made  for  those  trials. 


1.  INTRODUCTION 

The  current  research  project  addresses  the  need  of  the  US  Army  Medical  Research  and  Development 
Command  to  predict  dosages  and  deposition  rates  to  civilians  and  Indigenous  plants  and  wildlife 
resulting  from  the  use  of  military  smokes  In  field  training  exercises.  To  meet  this  need,  two  project 
objectives  are  set  forth:  (1)  carry  out  a  comprehensive  field  measurement  program  to  obtain  high  quality 
data  on  the  dosages  and  deposition  rates  resulting  from  smoke  emitted  by  specific  military  smoke 
generating  devices,  and  (2)  use  these  field  data,  combined  with  results  which  can  be  synthesized  from 
other  studie  of  snokes,  to  evaluate  the  performance  of  predictive  models.  Results  to  date  of  the  field 
effort  described  In  (I)  are  given  In  Ounn  et  al.  (1985). 
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The  Initial  effort  of  the  model  evaluation  program  has  focused  on  testing  model  predictions  with 
data  from  Smoke  Weeks  III  and  IV.  This  effort  not  only  provides  experience  with  the  use  of  the  models 
but  more  Importantly  highlights  critical  modeling  and  data  issues.  The  range  of  variability  among  the 
predictions  of  alternative  modeling  approaches  is  also  being  doci*ented.  Further,  sensitivity  tests 
will  Indicate  which  parameters  need  to  be  most  accurately  measured  In  our  field  studies.  The  effort 
will  establish  a  framework  for  the  evaluation  process  Itself  Including  the  use  of  graphical  and 
statistical  presentations  of  nodel/data  comparisons. 

This  paper  summarizes  the  early  efforts  under  this  model  evaluation  plan.  Consistent  with  overall 
project  objectives,  the  focus  has  been  on  the  evaluation  of  models  for  fog-oil  pltme  dispersion  In  flat 
terrain.  The  specific  data  sets  used  In  model  testing  are  Smoke  Meek  III,  Trials  4  and  19  along  with 
Smoke  Week  IV,  Trial  3.  Model/data  comparisons  were  made  for  concentration  dosages  along  one  (Smoke 
Week  IV)  or  two  (Smoke  Week  III)  horizontal  lines  of  sight  (LOS).  The  models  tested  were  COMBIC  (Hoock 
et  al.,  1982),  ACT  ll  (Sutherland  and  Hoock,  1982),  HAD  Puff  (Matlse  et  al.,  1982),  and  the  Ludwig  Model 
(Ludwig,  1977). 

Future  work  will  Involve  the  testing  of  other  U.S.  Army  models  Including  HECSMOKE-I  (Cheney  and 
Oumbauld,  1979),  HAZRD-2  (personal  communication,  July  1984,  Hr.  Ronald  Pennsyle,  Chemical  Research  and 
Development  Center,  Aberdeen  Proving  Ground,  Maryland),  SEMM  (Marchettl,  1979,  1980,  1981),  HoCaPD 
(Huang  and  Frost,  1982).  the  Ohmstede-Stenmark  (1981)  Model,  the  Ludwig  (1983)  Model,  and  MSMOKE 
(Hansen,  1984).  Models  developed  outside  the  U.S.  Army  community  which  show  special  promise  for  fog-oil 
dispersion  calculations  will  be  evaluated  as  well.  Models  with  such  unique  theoretical  features  are  the 
Ounn-Boughton  Model  (1984),  the  Rlso  Puff  Model  (Mlkkelsen  and  Larsen,  1984),  and  the  Petersen  (1984) 
Model.  All  models  will  be  tested  with  the  above  Smoke  Week  data  along  with  new  data  obtained  as  part  of 
the  ongoing  field  work. 


2.  MODEL  FORMULATIOHS 

In  this  paper,  four  models  are  tested  with  field  data  on  measured  concentration  dosages:  COMBIC, 
ACT  II,  MAD  Puff,  and  Ludwig  (1977).  ACT  II  and  MAD  Puff  were  provided  to  us  by  Optimetrics  Inc.  as 
part  of  the  ACTMAD  computer  model.  If  run  In  Its  ACT  II  mode.  It  is  identical  to  the  standard  ACT  II 
model  (Sutherland  and  Hoock,  1982)  with  only  a  few  minor  alterations.  If  run  In  Its  HAD  Puff  mode,  it 
uses  the  MAD  Puff  turbulent  transport  methodology  to  be  discussed  below.  The  CCMBIC  Model  was  obtaire: 
from  the  Atiospherlc  Sciences  laboratory  and  the  Ludwig  Model  was  obtained  from  Dr.  Ludwig  at  f 
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Stanford  Research  Institute.  The  Ludwig  (1977)  Model  Is  also  referred  to  as  the  Behavior  of 
Environmental  Aerosol  Releases  (BEAR)  Model  and  Is  the  most  appropriate  of  the  Ludwig  models  to  be 
tested  with  the  Smoke  Week  data. 

The  key  features  of  each  of  the  four  models  are  compared  in  Table  1.  As  may  be  seen,  the  COMBIC. 
ACT  II.  MAD  Puff,  and  Ludwig  Models  are  all  of  Gaussian  type.  In  COMBIC.  the  continuous  plume  (rather 
than  the  puff  formulation)  Is  applied  to  fog  oil  released  from  smoke  generators.  The  model  developers 
believe  that  releases  of  10  seconds  or  greater  are  best  simulated  with  a  continuous  plume  formulation. 
They  cite  the  fact  that  observed  fog-oil  plumes  do  not  have  a  leading  edge  or  trailing  portion  outside 
of  the  main  body  of  the  plume.  Except  for  the  lowest  wind  speeds,  longitudinal  advectlon  dominates 
longitudinal  diffusion.  As  a  result,  diffusion  In  the  longitudinal  direction  need  not  be  modeled.  The 
ACT  II.  MAD  Puff  and  Ludwig  Models  simulate  the  generator  effluent  as  a  series  of  discrete  puffs  during 
the  time  of  generator  operation.  The  advantage  to  a  puff  release  simulation  Is  that  a  time  varying 
release  can  be  more  accurately  treated. 


TABLE  1.  COMPARISON  OF  THEORETICAL  FEATURES  OF  COMBIC.  ACT  II.  MAD  PUFF. 
AND  LUDWIG  (1977)  MODELS  FOR  FOG-OIL  APPLICATIONS 


Model 

Approach 

Formulas  for  . 

V  v  and  °z 

Meteorol oglcal 
Data  Input 

Plume 

Rise 

Deposl  tlon 

Uses  Measured 
Turbul ence 
Quantl tl es . 

0^  and  3} 

COMBIC 

Gaussl  an 

PI  ume 

Continuous 

plume 

Steady  state 

,N0 

No 

No 

ACT  II 

Gaussl  an 

Puff 

Continuous 
pi  ume 

Steady  state 

No 

No 

No 

MAD  Puff 

Gaussian  Puff 
with  Random 
Trajectory 
Methodol  ogy 

Instantaneous 

puff 

Steady  state 

No 

No 

Yes 

Ludwl  a 
(1977) 

Gaussl  an 

Puff 

Continuous 
pi  ume 

Time  dependent 

Tes 

No 

No 

Botn  ACT  II.  COMBIC,  and  MAD  Puff  Models  reoulre  the  use  of  time-averaged  meteorological  data;  as  a 
result,  the  models  are  applied  to  a  pceudo  steady-state  meteorological  system.  For  all  three  models, 
the  time-averaging  Is  done  over  ?  period  of  operation  of  the  smoke  generator.  The  Ludwig  Model 
permits  as  input  time-dependent  rvteorologlcal  data  In  Its  puff  formulation.  Each  time  that  new 
meteorological  Inputs  are  read,  a  new  time  step  for  the  release  of  cuffs  Is  determined.  The  Ludwig 
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treatment  permits  the  simulation  of  changing  wind  speeds  and  directions  during  the  period  of  smoke 
generation.  As  a  result,  plume  meander  Is  treated  by  means  of  this  time-dependent  formulation. 

In  spite  of  the  fact  that  time- averaged  meteorological  data  are  used  In  HAD  Puff,  the  model  uses 
measured  turbulence  quantities  (o9  and  o#)  to  add  random  transport  components  to  the  trajectory  of  each 
puff.  In  the  HAD  Puff  mode,  the  puffs  are  moved  In  a  semi-random  fashion  based  upon  a  turbulent  fluid 
parcel  transport  theory  developed  by  S.R.  Hanna  (1979).  As  a  result,  each  puff  has  a  variable 
trajectory  based  on  the  simulation  of  atmospheric  turbulence.  The  growth  of  the  puffs  with  downwind 
distance  Is.  however,  based  on  the  plume  dispersion  parameters.  ox.  o^.  and  o2. 

C0MB1C  and  ACT  II  provide  smooth  Gaussian  plumes  with  a  straight-line  trajectory  for  the  continuous 
plume  (COMBIC)  and  puffs  (ACT  II).  MAD  Puff  and  Ludwig  create  Inhomogeneous  plumes  by  means  of  their 
variable  puff  trajectory  formulations.  MAD  Puff  uses  a  simulation  of  atmospheric  turbulence  to  provide 
the  randomness  whereas  the  Ludwig  Model  employs  time  dependent  wind  speeds  and  directions.  Figure  1 
provides  a  characterization  of  the  smoke  plumes  predicted  by  each  model  as  compared  to  a  "real"  smoke 
plume.  — . 


C0TB1C  ACT  II  HQ  Puff 


FIGURE  X.  SMOKE  PLUME  CHARACTERIZATION  OF  COMBIC,  ACT  II, 

MAD  PUFF,  AND  LUDWIG  (1977)  AS  COMPARED  TO  A 
TYPICAL  ■real'  PLUME 

Each  of  the  above  four  Gaussian  Models  uses  different  power-law  formulations  for  the  horizontal  and 
vertlcz!  dispersion  coefficients.  The  CCM3IC  Model  uses  the  a  and  oz  coefficients  recomenc^i  by 
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Hansen  (1979)  for  continuous  plumes.  In  this  notation,  x  is  downwind,  y  is  crosswind,  and  2  is 
vertical.  Oue  to  the  continuous  Gaussian  plume  assumption  for  fog-oil  plumes.  COMBIC  assumes  ax*Q.  For 
the  ACT  II  puff  release.,  the  continuous  plume  sigmas  from  the  Hansen  (1979)  reference  are  used  but  with 
ox«cy  Ludwig  uses  continuous  plume  sigmas  for  his  puffs;  is  assumed  with  variations  with 

downwind  distance  obtained  from  Turner's  workbook  (1969).  The  o2  values  used  by  Ludwig  are  a 
modification  of  Turner's  values  for  each  stability  class.  In  the  COMBIC  and  ACT  II  Models.  o2  is  a 
function  of  roughness  length  in  addition  to  stability  class. 

The  plume  sigmas  used  by  the  COMBIC.  ACT  II.  and  Ludwig  Models  were  obtained  fra*  data  on  long-term 
time  averages  of  steady-state  plumes;  these  sigmas  include  the  effects  of  the  low  frequency  wandering  of 
the  plume  as  well  as  small-scale  diffusion.  The  MAD  Puff  Model  attempts  to  simulate  the  meandering  of 
the  plume  separately  from  small-scale  diffusion  processes.  As  a  result,  the  sigmas  used  by  MAO  Puff 
represent  those  of  Instantaneous  puffs  where  essentially  the  effects  of  the  low  frequency  meandering 
have  been  removed  from  the  sigmas  obtained  from  long-term  averaging.  The  formulas  for  the 
Instantaneous  puff  sigmas  are  also  obtained  fra»  Hansen  (1979).  From  Hansen,  the  growth  rate  of  the 
Instantaneous  puff  sigmas  is  simply  2/3  the  growth  rate  of  the  long-term  sigmas  for  the  y  and  2 
directions  with  a  factor  of  0.74  used  for  the  a  growth  formula.  As  a  result,  o  Is  greater  than  0  In 
the  MAO  Puff  formulation.  Clearly,  the  different  simulations  of  dispersion  among  the  models  will  lead 
to  differences  in  model  predictions. 

The  COMBIC.  ACT  II  and  Ludwig  models  assume  complete  reflection  of  the  plume  at  the  ground 
surface.  The  option  for  complete  reflection  was  chosen  In  the  fog-oil  runs  for  MAD  Puff.  Mo  deposition 
of  fog  oil  Is  also  assumed.  The  Ludwig  Model  has  a  built-in  treatment  for  plume  rise  wtiereas  the  user 
In  COMBIC.  ACT  II.  and  MAD  Puff  suppresses  the  treatment  of  plume  rise  for  fog-oil  applications  with  his 
choice  of  Input  parameters. 


3.  SMOKE  WEEK  III  AMD  IV  DATA 


Three  trials  were  chosen  for  the  purpose  of  comparing  model  predictions  with  data.  These  are  Smoke 
Week  III.  Trials  4  and  19.  and  Smoke  Week  IV.  Trial  3.  These  trials  employed  the  M3A3  smoke 
generator.  Data  from  these  trials  were  obtained  largely  from  the  Atmospheric  Aerosol  ana  Optics  Data 
Library  ( AA00L )  of  Science  and  Tech.ology  Corporation  with  supporting  Information  from  0PM  Smk/Cbs 
(personal  communication.  August  19B4,  Mr.  Howard  Smalley,  Aberdeen  Proving  Ground,  Maryland).  Burgess 
and  Nielsen  (1  J ) .  and  Nelson  et  al .  (1982). 
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LOS  2:  33  APs  (Approx.  9  m  apart) 

101  CIs  (Approx.  3  h  apart) 
LOS  3:  9  APs  (Approx.  18  m  apart) 

27  CIs  (Approx.  6  m  apart) 
Instrumentation  layout  centered 
about  met  toner  (grid  origin) 


WIND  DIRECTION 
(206  DEG.  E.  OF  N.) 


32 -METER  MET  TOWER 


SOURCE 


SCALE:  "ETFBS 


50  100 


FIGURE  2.  CONFIGURATION  OF  SOURCE.  RECEPTORS.  AND  METEOROLOGICAL 
TOWER  FOR  SMOKE  WEEK  III.  TRIAL  4 


LOS  2:  33  APs  (Approx.  9  h  apart)  j 

101  CIs  (Approx.  3  m  apart) 
LOS  3:  9  APs  (Approx.  18  m  apart) 

27  CIs  (Approx.  6  m  apart) 
Instrumentation  layout  centered 

ABOUT  MET  TOMER  (GRID  ORIGIN) 


WIND  DIRECTION 
(196  DFG.  F.  OF  N.) 


32-METER  MET  TOWER 


SOURCE 


SCALE:  METERS 


FIGURE  3.  CONFIGURATION  OF  SOURCE.  RECEPTORS.  AND  METEOROLOGICAL 
TOWER  FOR  SMOKE  WEEK  III.  TRIAL  19 
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LOS  1:  5  Pairs  of  C!s  (Approx.  18  m  apart) 
33  APs  (Approx.  9  m  apart) 
Instrumentation  layout  centered 

ABOUT  HET  TOWER  (GRID  ORIGIN) 


wind  niRErrioi 
(119  DEG .  E.  OF  N.) 
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FIGURE  4.  CONFIGURATION  OF  SOURCE.  RECEPTORS.  AND  METEOROLOGICAL 
TOWER  FOR  SMOKE  WEEK  IV.  TRIAL  3 


In  addition,  four  2-m  meteorological  masts  located  along  the  main  sampling  line  were  equipped  to  measure 
horizontal  wind  direction  and  wind  speed.  Each  of  the  models  used  only  wind  speed  and  direction  at  a 
10-m  elevation  obtained  by  interpolating  the  8  and  16-m  wind  data  on  the  meteorological  tower.  The  wind 
profile  exponent  computed  from  the  measured  wind  profile  was  used  as  an  Input  parameter  by  the  ACTMAD 
Model.  PasQuill  stability  categories  were  determined  using  the  Turner  metnod  based  on  wind  speed  and 
solar  radiation  Index,  modified  by  cloud  cover  and  cloud  height.  None  of  the  mast  data  or  Individual 
wind  measurements  at  different  tower  levels  could  be  Input  to  the  models. 

The  Smoke  Week  IV  tests  were  conducted  at  Test  Area  1  ( TA- 1 )  at  Redstone  Arsenal.  Huntsville. 
Alabama.  TA-1  Is  a  relatively  flat,  open  field  approximately  6  km  long  and  1-2  km  wide.  Part  of  the 
test  area  was  damp,  and  the  lowest  areas  contained  pools  of  water.  TA-1  Is  surrounded  by  dense  hardwood 
forest.  Tree  heights  are  typically  10  m.  In  the  Smoke  Week  IV  tests,  only  one  line  of  samplers  was 
used  associated  with  LOS  1  of  the  total  of  four  optical  lines  of  sight.  Figure  4  provides  a  sketch  of 
this  sampling  line,  the  45-m  meteorological  tower,  and  the  generator  location.  On  the  sampling  line 
were  33  aerosol  photometers  and  12  aerosol  samplers  (chemical  Implngers  and/or  filter  samplers). 
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Samplers  were  positioned  In  pairs  at  five  points  along  LOS  1  with  two  on  a  mobile  measurement  unit, 

1.5  m  above  ground.  The  sampling  line  also  had  33  aerosol  photometers  at  9-m  Intervals. 

A  45-a  meteorological  tower  was  located  at  the  center  of  the  sampling  line  and  was  equipped  to 
measure  the  following: 

(a)  temperature  and  dew  point  at  the  2  and  32-m  levels, 

(b)  temperature  at  the  0.5,  2,  4,  10,  16,  32,  and  45-m  levels, 

(c)  horizontal  components  of  the  wind  direction  and  speed  at  the  2,  4,  10,  16,  32,  and  45-m  levels, 
and 

(d)  vertical  component  of  wind  direction  at  the  10  and  32-«  levels. 

o0  was  available  at  the  2.  4,  10,  16,  32,  and  45-a  levels  while  was  available  only  at  the  10  a-d  12-m 
levels.  A  wind  profile  exponent  required  by  ACT  HAD  was  available  from  the  wind  data.  Atmospheric 
Stability  was  provided  In  terms  of  Pasqulll  categories  determined  as  before  from  sky  conditions  and  the 
10-m  wind  speed. 


4.  PREPARATION  OF  MODEL  INPUTS 

As  noted  above,  each  model  has  special  Input  requirements  which  utilize  a  subset  of  the  available 
data.  Considering  that  COMBIC,  ACT  II,  and  MAD  Puff  require  meteorological  data  that  are  averaged  over 
the  generator  run  time,  their  Input  data  are  quite  similar.  Table  2  lists  the  key  Input  to  those 

models.  All  models  Including  that  of  Ludwig  require  wind  speed  and  direction  at  the  10-m  level.  The 

Ludwig  Model  Is  different  In  that  It  requires  a  time-dependent  Input.  Meteorological  data  were  Input  to 
that  model  In  five  second  Intervals.  To  be  noted  from  Table  2  Is  that  all  three  cases  happen  to 

repretent  stability  class  C.  Also,  an  assunptlon  of  40  gal/hr  of  fog  oil  release  from  the  generator  had 
to  be  assumed  for  Smoke  Week  III,  Trials  4  and  19.  No  data  on  fog  oil  use  were  reasured  In  those 

experiments.  The  data  from  Smoke  Week  IV,  Trial  3  reveal  only  a  24  gal/hr  release  rate.  Generator 
operation  times  were  short  varying  from  only  5  to  11  minutes.  Longer  periods  of  release  are  of  Interest 
for  health  effects  research;  such  longer  period,  would  likely  have  the  wind  vary  In  direction  much  more 
and  lead  to  a  greater  challenge  for  the  models. 
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TABLE  2.  SUMMARY  OF  KEY  MODEL  INPUTS  FOR 
LUDWIG  (1977)  MODELS  AS  APPLIED 
OATA  SETS 

C0MBIC,  ACT  II,  MAD  PUFF  AND 

TO  THE  SMOKE  WEEK  III  AND  IV 

Input 

Parameter 

Smoke  Week  III 
Trial  4 

Smoke  Week  III 
Trial  19 

Smoke  Week  IV 
Trial  3 

Wind  speed  (10  m),  m/s 

5.0 

3.7 

4.0 

Wind  direction  (10  m),  deg 

206 

196 

119 

Wind  speed  power  law 
exponent 

0.11 

0.14 

0.116 

Stability  class 

C 

C 

C 

Roughness  height,  cm 

50 

50 

10 

Generator  release  rate, 
gal /hr 

Not 

measured* 

Not 

measured* 

23.8 

Generator  run  time,  rtrin 

7 

11 

5 

Location  of  X  axis 
(line  of  samplers), 
deg  east  of  north 

109 

108 

74 

Grid  origin  (0,0) 

32 -m  met  tower 

32-m  met  tower 

45-m  met  towe  - 

Generator  coordinates,  m 

(-30,  -60) 

(-15.  -50) 

(90,  -80) 

*It  was  assumed  that  the  nominal  value  of  40 
fog  ol  1 . 

gal/hr  was  the  release  rate  of 

Special  Parameters  (for  MAD 

Puff)  only: 

Friction  velocity, 
m/s 

0.7 

0.7 

0.5 

oe,  deg  (at  4-m  level ) 

12.6 

13.6 

8.0 

n?,  deg  (at  4-m  level ) 

6.0 

8.2 

4.5 

(at  10-m  level 

computer  code  was  modified 

to  permit  i_he 

calculation  of 

dosages  from 

concentration  variation  with  time.  C0M8IC,  ACT  II,  ana  MAD  Puff  were  originally  developed 
obscuration  calculations.  The  transport  and  diffusion  portions  of  the  code  are  relatively 
permit  the  computation  of  concentrations  with  time  at  any  receptor  from  which  dosages 
concentrations  with  time)  were  calculated. 


p  re  a ' 

to  ; r; 
■  small 
( 1  n  t  e  a  ’ 
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All  computer  runs  were  carried  out  on  a  Ridge  32  mini -computer.  There  Is  roughly  a  factor  of  6 
greater  computation  time  for  the  Ridge  as  compared  to  an  IBM  3033  mainframe.  COHBIC  was  run  with  one 
second  time  steps  and  ACT  TI  and  MAO  Puff  were  run  with  one  second  between  puff  releases.  The 
equivalent  run  times  of  the  COHBIC,  ACT  II.  and  MAD  Puff  models  for  an  IBM  3033  were  about  10  minutes 
for  each  line  of  sight.  Each  of  the  models  was  run  within  Its  original  computer  code  which  carried  out 
many  other  calculations  (such  as  transmittance  across  lines  of  sight)  that  were  largely  irrelevant  to 
the  objectives  of  the  present  validation  exercise.  A*-  a  result,  the  run  times  are  considerable  larger 
than  If  the  models  were  ccdca  solely  for  transport  ana  diffusion  calculations  with  the  subsequent  output 
of  dosages.  The  Ludwig  Model  was  run  with  10  seconds  between  puff  releases;  the  computer  code  ran  an 
equivalent  of  only  30  seconds  on  the  IBM  3033.  This  model  was  set  up  directly  to  produce  concentration 
predictions  with  time  from  which  we  prepared  dosage  values. 

5.  RESULTS  OF  M00EL/DATA  COMPARISONS 

The  comparison  of  dosage  predictions  and  data  for  the  three  Smoke  Week  trials  are  given  In 
Figures  5-9.  Smoke  Week  III.  Trial  4  comparisons  for  LOS  3  and  2  are  given  respectively  In  Figures  5 
and  6.  For  Smoke  Week  III.  Trial  19.  comparisons  for  dosage  are  presented  for  LOS  3  and  2  In  Figures  7 
and  8.  LOS  1  predictions  and  data  for  Smoke  Week  IV.  Trial  3  are  given  In  Figure  9. 

Examination  of  the  comparisons  reveals  the  following  observations: 

(a)  For  the  short  distances  (approximately  50-150  m)  Involved,  the  predictions  of  the  models  are 
very  Similar  to  each  other  and  are  within  factors  of  2-3  of  the  data.  This  flnd'ng  Is  not  surprising 
considering  the  simlla'  theoretical  formulations  of  the  models.  It  Is  not  clear  whether  the  similarity 
In  predictions  remains  at  larger  distances.  This  latter  question  will  be  Investigated. 

(b)  All  model  predictions  of  the  peak  dosage  are  offset  from  the  data.  This  result  Is  to  be 
expected  for  CCM8IC  and  ACT  II  since  these  models  employ  time-averaged  meteorological  data  and.  thereby, 
cannot  resolve  wind  direction  variations  during  the  period  of  transport.  The  MAD  Puff  Model  accounts 
for  wind  direction  fluctuations  In  Its  u'e  of  o?  and  a  and  the  Ludwig  model  handles  wind  direction 
variations  through  Its  Input  of  time-dependent  meteorology.  The  Ludwig  Model  performs  best  In  te~ns  of 
estimating  the  location  of  the  peax  dosage. 
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SMOKE  WEEK  IV  -  TRIAL  3  —  LOS1  DOSAGE 


FIGURE  9.  COMPARISON  OF  DOSAGE  PREDICTIONS  OF  COMBIC.  ACT  II.  MAD  PUFF. 
ANO  LUDWIG  (1977)  WITH  FIELD  DATA  FROM  SMOKE  WEEK  IV. 

TRIAL  3.  LOS1 


(c)  For  short  release  tines  and  the  calculation  of  dosages,  the  randomization  of  puff  positions  (as 
In  MAD  Puff)  do  not  appear  to  Improve  model  predictions.  Comparisons  for  further  distance,  will  be  made 
In  future  tests. 

(d)  In  general,  the  models  tend  to  underestimate  lateral  spreading.  This  disagreement  with  the 
data  nay  be  the  result  of  Inadequate  treatment  of  wind  direction  fluctuations  during  the  period  of  smoke 
release,  and 

(e)  The  lack  of  data  on  the  fog-oil  release  rate  for  the  Smoke  Week  III  tests  created  some 
uncertainty  In  model  Inputs.  Apparently,  the  nominal  40  gal/hr  rate  used  appeared  to  be  a  reasonable 
value;  however,  more  enact  comparisons  of  models  to  data  require  measurements  of  the  rate  of  release  of 
fog  oil  wlti  time  as  well  as  the  total  fog  oil  consumed  during  generator  run  time.  The  models  performed 
best  with  the  data  of  Smoke  Week  IV,  Trial  3;  for  this  trial,  the  consumption  of  fog  oil  during  the  test 
was  measured. 
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The  data  from  the  Smoke  Week  trials  are  very  limited  In  quantity  and  distance  of  coverage.  In  fact,  the 
three  trials  used  In  this  paper  we-e  the  best  available  from  the  Smoke  Week  tests  yet  they  happened  tc 
have  very  similar  meteorological  conditions  (stability  class  C.  wind  speeds  of  about  4  m/s).  Additional 
data  are  required  In  order  to  make  more  definitive  conclusions  on  model  performance. 

Future  work  will  Include: 

(1)  the  testing  of  new  models  with  Smone  Week  data.  Additional  Smoke  Week  trials  using  fog  oil 
will  be  used;  however,  most  of  the  other  trials  Involve  the  movement  of  the  smoke  generator  during  the 
test.  Some  models  cannot  treat  this  change  in  source  location  effectively  within  their  theoretical 
framework. 


(II)  an  evaluation  of  the  magnitude  of  differences  In  model  predictions  at  distances  In  the  range 
500-2000  m  from  the  smoke  generator.  It  Is  to  be  determined  whether  there  are  greater  differences  In 
model  predictions  than  were  shown  here  at  about  50-150  m. 

(III)  a  sensitivity  analysis  of  the  models  to  key  Input  parameters  to  determine  which  ones  require 
the  greatest  attention  In  the  associated  field  program. 

( 1 v )  the  development  of  a  ’generic"  model  In  which  numerous  different  alternative  approaches  can  be 
tested  within  one  code.  In  this  way.  the  fundamental  differences  between  modeling  approaches  can  be 
compared  Irrespective  of  the  details  of  the  specific  model  formulation,  and 

(v)  the  testing  of  all  models  cited  above  with  new  fog  oil  trials  presented  In  Dunn  et  al . 
(1985). 
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